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1.0   INTRODUCTION 

Since  the  development  of  the  scintillation  detector  and  the  pulse- 
height  analyzer  a  significant  amount  of  work  has  been  done  in  the  area  of 
gamma-ray  spectroscopy.   However,  the  pulse-height  spectrum  obtained  through 
the  use  of  such  a  system  is  not  a  proportional  representation  of  the  energy 
spectrum  of  gamma  rays  striking  the  scintillation  detector.   This  is  the 
result  of  distortions  in  both  the  scintillation  detector  crystal  and  the 
associated  electronic  equipment,  together  comprising  the  gamma-ray  spectrom- 
eter system. 

With  the  realization  that  the  experimentally  observed  spectrum  is  not 
an  accurate  representation  of  the  true  energy  spectrum,  a  number  of  mathe- 
matical techniques  (1-9)  have  been  proposed  to  enable  the  determination  of 
the  true  gamma-ray  energy  spectrum  from  the  experimental  spectrum.  These 
techniques  are  referred  to  as  unfolding  and  since  they  usually  require  the 
use  of  high  speed  computers  they  are  commonly  known  as  unfolding  codes. 

Although  each  of  these  unfolding  techniques  has  its  merit  there  are 
variations  in  the  mathematical  approaches  used  in  arriving  at  the  solutions 
as  well  as  the  physical  models  on  which  they  are  based.   Hence  the  applica- 
tion of  a  particular  unfolding  technique  may  vary  with  experimenters  because 
of  dissimilarities   in  scintillation  spectrometer  systems.   Therefore,  even 
though  a  general  unfolding  technique  is  used,  each  experimenter  must  develop 
an  unfolding  code  which  will  accurately  describe  the  distortions  introduced 
into  the  experimentally  observed  spectrum  by  the  scintillation  spectrometer 
system  and  also  take  into  account  the  experimental  environment  in  which  the 
measurements  were  made. 


Many  papers  (1,2,4,6,7,8)  have  been  published  listing  the  results  of 
the  applications  of  unfolding  codes  for  a  variety  of  detection  systems  and 
experimental  environments.   However,  there  has  not  been  available  a  generally 
accepted  method  for  determining  how  well  the  unfolded  spectrum  represents 
the  true  energy  spectrum  of  radiation  incident  upon  the  detector.   Unfolding 
codes  are  often  confirmed  by  applying  them  to  a  "standard  complex  spectrum." 
A  spectrum  of  this  type  is  generated  by  combining  spectra  resulting  from 
standard  monoenergetic  sources.   Then  the  gamma-ray  energies  and  their 
intensities  which  are  obtained  by  unfolding  the  "standard  complex  spectrum" 
can  be  compared  to  the  energies  and  intensities  of  the  standard  monoener- 
getic sources.   The  problem  with  using  such  a  method  is  that  the  unfolding 
codes  are  based  on  spectrometer  system  responses  to  standard  monoenergetic 
sources  such  as  used  to  generate  the  "standard  complex  spectrum."   There- 
fore this  procedure  does  not  really  test  the  validity  of  the  unfolding  code 
for  gamm.a  rays  whose  energies  are  different  from  those  of  the  standard 
monoenergetic  sources.   In  practice  such  spectra  are  often  encountered.   In 
addition,  the  energies  of  the  gamma  rays  incident  upon  the  detector  are 
often  nearly  equal  and  therefore  their  individual  contributions  to  the  com- 
plex spectrum  can  not  readily  be  determined.   In  this  case  one  does  not 
know  either  the  number  or  the  energies  of  the  gamma  rays  which  contribute 
to  the  complex  spectrum  and  therefore  it  is  difficult  to  theoretically 
predict  the  shape  and  intensity  of  the  true  energy  distribution. 

The  purpose  of  this  work  has  been  to  determine  experimentally  the 
validity  of  using  an  unfolding  code,  developed  for  a  collimator-spectrometer 
system,  as  a  tool  in  the  computation  of  exposure  doses  resulting  from 
experimentally  measured  complex  spectra.   This  determination  is  based  on  a 


comparison  of  exposure  doses  computed  on  the  basis  of  measured  complex 
gamma-ray  energy  spectra  with  exposure  doses  measured  with  an  air  equiva- 
lent ionization  chamber.   Although  this  report  is  concerned  with  the 
accuracy  of  exposure  doses  calculated  using  an  unfolding  code,  a  detailed 
analysis  of  the  unfolding  code  in  question  is  not  intended. 


2.0   THEORETICAL  DEVELOPMENT 
2.1  Unfolding  Technique 

In  the  design  of  an  experiment  to  determine  the  accuracy  and  thereby 
obtain  a  measure  of  the  validity  of  an  unfolding  technique,  the  first  re- 
quirement is  a  basic  understanding  of  the  technique  in  question.  Hence  it 
is  imperative  to  know  under  what  experimental  conditions  the  technique  is 
applicable,  the  general  principles  involved  and  the  form  in  which  results 
are  obtained  as  well  as  what  they  represent. 

The  unfolding  technique  under  investigation  in  this  report  was  developed 
by  J.  A.  Baran  (10)  as  a  result  of  a  research  project  sponsored  by  the 
Department  of  Defense  Office  of  Civil  Defense  under  Contract  OCD-OS-63-74. 
The  primary  purpose  of  this  technique  was  to  provide  a  method  for  determining 
the  true  gamma-ray  energy  spectrum  from  an  experimentally  measured  pulse- 
height  distribution.   A  knowledge  of  the  true  gamma-ray  energy  spectrum  was 
necessary  for  the  computation  of  differential  exposure  doses  resulting  from 
the  scattering  of  gamma  rays  from  small  areas  on  a  concrete  slab.   Scattering 
of  gamma  rays  from  concrete  results  in  an  emergent  flux  of  gamma  rays  with  a 
distribution  of  energies,  i.e.,  a  plane  source  of  poliJ-energetic  gamma  rays. 

The  experimental  apparatus  for  which  this  technique  was  developed  con- 
sisted of  a  scintillation  detector  located  within  a  collimator  and  a  pulse- 
height  analyzer.  Hereinafter  this  group  of  Instruments  will  be  referred  to 
as  the  collimator-spectrometer  system.  Since  the  source  to  which  the 
detector  was  exposed  consisted  of  scattered  gamma  rays,  the  unfolding  tech- 
nique was  formulated  to  analyze  complex  gamma-ray  spectra  over  a  wide  range 
of  gamma-ray  energies. 


A  detailed  development  of  this  unfolding  technique  is  beyond  the  scope 
of  this  report;  however,  a  complete  analysis,  including  computer  codes,  is 
available  in  reference  10.   For  the  purposes  of  this  report,  the  problem  of 
unfolding  can  be  simply  stated  through  the  use  of  the  matrix  equation 

C  =  R  N  (1) 

where     C   is  the  experimentally  observed  spectrum, 
R  is  a  response  matrix  and 
N   is  the  desired  unfolded  pulse-height  distribution. 

In  the  above  equation,  C  and  N  are  column  vectors  and  R  is  a  square  matrix. 
The  response  matrix,  which  relates  the  experimental  spectrum  to  the  unfolded 
pulse-height  distribution,  corrects  the  experimental  spectrum  for  the  dis- 
tortions introduced  by  the  collimator-spectrometer  system,  i.e.,  the  scin- 
tillation detector,  the  collimator  and  the  associated  electronics  used  in 
making  the  experimental  measurements.   This  matrix  was  obtained  by  consider- 
ing the  response  of  the  system  to  a  number  of  calibrated  monoenergetic  gamma- 
ray  sources  over  the  energy  range  of  interest. 

Solving  Eq.  (1)  for  the  unfolded  pulse-height  distribution  results  in 

N  =  R-1  C  .  (2) 

The  unfolded  pulse-height  distribution  obtained  from  this  equation  is  thus 
a  measure  of  the  energy  distribution  of  the  gamma-ray  flux  at  the  center  of 
the  crystal.   It  should  be  noted  that  this  pulse-height  distribution  has 
been  corrected  for  the  effects  of  the  collimator,  i.e.,  the  distortion  in 
the  collimator,  as  well  as  the  response  of  the  spectrometer  system. 


2.2  Design  of  an  Appropriate  Experiment 

With  the  preceding  discussion  concerning  the  conditions  under  which  the 
unfolding  technique  is  applicable  and  an  understanding  of  what  the  results 
represent,  an  appropriate  experiment  can  be  designed  to  determine  the 
accuracy  of  this  technique.   However,  at  this  point,  it  is  advantageous  to 
examine  this  problem  in  general  before  considering  the  design  of  an  experi- 
ment for  this  particular  technique. 

In  general  the  problem  of  determining  the  accuracy  of  a  particular 
technique  which  is  to  be  used  in  the  calculation  of  a  specific  experimentally 
observable  physical  phenomenon  often  results  in  a  pseudo  calibration  of  the 
technique  in  question.   When  confronted  with  such  a  problem,  a  logical  solu- 
tion is  to  design  an  experiment  such  that  the  quantity  to  be  measured  is 
either  known  precisely  or  can  be  obtained  directly  using  another  method  whose 
accuracy  has  previously  been  established.   Assuming  that  such  an  experiment 
can  be  designed  and  carried  out,  then  a  comparison  of  the  value  which  was 
determined  experimentally,  using  the  technique  under  investigation,  to  the 
known  or  accepted  value  will  be  a  measure  of  the  accuracy  of  the  procedure 
in  question. 

However,  it  is  often  found  that  in  the  use  of  such  an  experiment  to 
determine  the  accuracy  of  a  specific  technique,  the  value  of  the  quantity 
to  be  measured  experimentally  neither  is  known  precisely  nor  can  it  be 
calculated  directly  from  theoretical  considerations  alone.   \\Tien  this  is 
the  case  a  generally  accepted  course  of  action  is  to  determine  the  standard 
or  accepted  value  by  experimentally  measuring  the  same  quantity  using  a 
completely  independent  method,  whose  accuracy  can  be  well  established. 


In  order  to  determine  the  accuracy  of  the  unfolding  technique,  the 
experiment  to  be  designed  must  be  such  that  the  resulting  confidence  on  this 
technique  will  be  valid  under  the  conditions  for  which  it  is  to  be  used. 
Thus  a  complex  gamma-ray  spectrum  must  be  investigated  if  this  technique  is 
to  be  properly  evaluated.   Any  theoretical  evaluation  of  the  unfolding 
technique  would  therefore  require  a  knowledge  of  the  gamma-ray  energy  dis- 
tribution under  investigation.   Since  calibrated  plane  sources  emitting 
gamma  rays  over  a  wide  energy  range  are  nonexistent,  such  plane  sources  must 
necessarily  be  simulated.   The  simulation  of  such  a  source  would  then  make  it 
necessary  to  calculate  the  energy  distribution  of  the  gamma-ray  flux  under 
investigation  and  such  a  calculation  would  be  quite  involved,  and  then  only 
an  approximation.   Thus  it  is  advantageous  to  determine  the  standard  value 
experimentally.   These  values  could  then  be  compared  to  the  results  of  the 
collimator-spectrometer  measurements  in  order  to  determine  the  accuracy  of 
the  unfolding  technique. 

Hence  an  experiment  designed  to  fulfill  the  objectives  of  this  report 
must  meet  the  following  requirements: 

1.  a  plane  source  emitting  gamma  rays  over  a  spectrum  of  energies 
must  be  used  in  the  experimental  measurements, 

2.  measurements  must  be  made  using  several  different  plane  sources, 

3.  a  detector  must  be  selected  vjhich  can  be  used  in  the  same  experi- 
mental configuration  as  the  scintillation  detector,  and 

A.   the  quantities  measured  using  the  detector  of  3.  must  be  comparable 
to  those  obtained  from  the  unfolding  technique. 
Requirement  1.  is  necessary  in  order  that  a  nontrivial  determination  of  the 
accuracy  is  made  and  the  second  is  primarily  for  the  purpose  of  revealing  any 


systematic  errors.   The  detector  selected  to  be  used  in  the  independent 
measurements  must  be  of  reputable  accuracy,  must  have  physical  dimensions 
compatible  with  the  collimator  and  must  result  in  a  quantity  which  can  be 
directly  compared  to  the  results  obtained  using  the  unfolding  technique. 

An  air  equivalent  ionization  chamber  (dosimeter)  which  can  be  accu- 
rately calibrated,  obtained  in  a  variety  of  physical  sizes,  and  used  to 
measure  a  quantity  directly  proportional  to  the  exposure  dose,  meets  all  of 
the  previously  mentioned  requirements  for  the  detector  to  be  used  in  the 
measurements.   In  addition,  the  response  of  the  dosimeter  is  linear  with 
respect  to  the  dose  received.   Although  the  quantity  measured  by  the  dosimeter 
is  proportional  to  exposure  dose  and  that  obtained  from  the  spectral  measure- 
ments is  a  pulse-height  distribution,  both  measurements  may  readily  be  con- 
verted to  exposure  dose.   In  fact,  the  prime  objective  of  this  work  was  to 
determine  how  accurately  exposure  doses  could  be  determined  using  the  unfold- 
ing technique  to  analyze  the  spectral  data  obtained  using  the  collimator- 
spectrometer  system. 

Thus  the  experiment  used  to  determine  the  accuracy  of  the  unfolding 
technique  will  involve  the  comparison  of  exposure  doses  obtained  from  dosim- 
eter and  collimator-spectrometer  measurements.   Both  methods  of  measuring 
the  exposure  dose  must  be  carried  out  using  the  same  experimental  geometry 
and  plane  source.   The  dosimeter  measurements  must  necessarily  be  carried 
out  with  care  in  order  that  the  exposure  dose  obtained  can  be  accepted  with 
confidence  as  a  correct  measure  of  the  exposure  dose. 


2.3   Calibration  of  the  Dosimeter 

Since  the  measure  of  the  accuracy  of  the  unfolding  technique  can  be 
only  as  good  as  the  dosimeter  measurements,  the  dosimeter  must  be  accurately 
calibrated.   As  the  dosimeter  is  an  indirect  reading  instrument,  calibration 
of  the  dosimeter  is  necessary  if  an  exposure  dose  is  to  be  obtained. 

In  the  preceding  section  it  was  stated  that  the  dosimeter  and  collimator- 
spectrometer  measurements  must  be  made  using  the  same  plane  sources  and 
experimental  configurations.   This  implies  that  the  dosimeter  must  be  cali- 
brated while  located  within  the  collimator  and  exposed  to  a  plane  source  of 
gamma  radiation. 

Thus  the  calibration  curve  must  relate  the  dosimeter  response,  which  is 
in  microamperes  (ya) ,  to  an  exposure  dose  in  milliroentgens  (mr)  based  on  an 
exposure  to  a  plane  source  emitting  gamma  rays  with  an  energy  distribution 
identical  to  that  which  will  be  encountered  in  the  experiment  to  determine 
the  accuracy  of  the  unfolding  technique.   In  addition  the  calibration  must 
take  into  account  the  effects  of  the  collimator  so  that  the  dosimeter  results 
can  be  directly  compared  to  the  exposure  doses  obtained  from  the  collimator- 
spectrometer  measurements. 

Such  a  calibration  would  normally  be  performed  by  exposing  the  dosimeter 
to  a  calibrated  source  of  radiation  and  then  plotting  the  measured  dosimeter 
response  versus  the  theoretically  computed  exposure  dose  which  would  be 
expected  at  the  center  of  the  dosimeter.   However,  since  calibrated  plane 
polyenergetic  sources  are  not  readily  available,  as  was  previously  mentioned, 
a  method  of  calibrating  the  dosimeter  using  a  calibrated  monoenergetic  plane 
source  is  necessary.   This  calibration  must  incorporate  within  it  corrections 
for  the  effects  of  the  collimator  and  the  energy  dependence  of  the  dosimeter. 
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First  assume  that  it  is  possible  to  obtain  a  calibration  of  the  dosimeter 
in  air  and  exposed  to  a  calibrated  monoenergetic  point  source  which  emits 
gamma  rays  with  energy  Ei.   Thus  a  relationship  of  the  form 

VRlhl^a)   =  fe„  Vlhrnv)  (3) 

m         L]^   in 

El 
can  be  obtained,  where  PR   (ya)  is  the  experimentally  measured  dosimeter 

m  -^ 

El 
response,  P   (mr)  is  the  theoretical  exposure  dose  which  would  be  expected 

at  the  center  of  the  dosimeter  and  k        is  the  proportionality  constant  re- 
lating the  preceding  two  quantities. 

Now  let  a  -^  be  an  angular  correction  factor  defined  as  the  ratio  of 

o 

the  dosimeter  response  for  radiation  of  energy  E^  incident  normal  to  the  axis 
of  the  dosimeter  to  its  response  for  radiation  of  energy  E^  incident  at  some 

angle  9  to  the  axis  of  the  dosimeter.   Further  assume  that  the  dosimeter 

El 
calibration  resulting  in  Eq.  (3)  was  conducted  such  that  a  -"^  is  unity. 

The  dosimeter  must  next  be  placed  in  the  collimator  and  exposed  to  a 
calibrated  plane  source  of  radiation.   Ideally  this  plane  source  would  emit 
gamma  rays  with  an  energy  distribution  identical  to  that  which  \\rould  be 
investigated  in  the  collimator-spectrometer  measurements;  however,  since 
calibrating  the  dosimeter  in  the  collimator  requires  a  theoretical  calcula- 
tion of  the  exposure  dose,  this  is  not  feasible.   Also  such  a  calibrated 
plane  source  is  not  available.   Therefore  a  calibrated  monoenergetic  plane 
source  must  be  used  to  calibrate  the  dosimeter  in  this  geometric  configura- 
tion and  appropriate  corrections  made  to  account  for  the  fact  that  the  calibra- 
tion curve  is  to  be  used  in  the  measurement  of  radiation  emitted  from  a 
polyenergetic  plane  source  rather  than  a  monoenergetic  plane  source. 
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Consider  the  following  terms  which  will  be  used  in  describing  the 

calibration  of  the  dosimeter  located  in  the  collimator  and  exposed  to  a 

monoenergetic  plane  source: 

El 
-DR  ^  (ya)  =  experimentally  measured  dosimeter  response  to  a 

6   m 

monoenergetic  plane  source  emitting  gamma  rays  of 
energy  £■[  which  are  incident  on  the  dosimeter  at  an 
angle  6  to  the  axis  of  the  dosimeter, 
D  J^(mr)  =  theoretical  exposure  dose  resulting  from  exposure 

to  the  plane  source  emitting  gamm.a  rays  of  energy  E^ . 

V 

The  quantity  ^DR  '■  (ya)  is  the  response  of  the  dosimeter  due  to  direct  radia- 

^      -^  6  m  "^ 

ticn  plus  that  radiation  which  reaches  the  detector  as  a  result  of  inter- 

Ei 
actions  with  the  collimator.   However,  D  Mmr)  represents  the  exposure  dose 

which  would  be  expected  as  a  result  of  direct  radiation  alone.   Direct 

radiation  is  considered  to  be  those  gamma  rays  which  are  emitted  from  the 

source  in  such  a  direction  that  they  can  reach  the  detector  without  the  need 

of  an  intermediate  scatter.   Therefore,  when  discussing  the  plane  source 

measurements  using  the  collimator,  direct  radiation  refers  to  that  radiation 

emanating  from  the  disk  on  the  source  plane  subtended  by  the  solid  angle  of 

El 
the  collimator.   The  theoretical  dosimeter  response  D  Mya)  resulting  from 

o  ih 

exposure  of  the  dosimeter,  located  within  the  collimator,  to  a  monoenergetic 
plane  source  emitting  gamma  rays  of  energy  E^  which  are  incident  on  the 
dosimeter  at  an  angle  6  to  the  dosimeter  axis,  can  be  computed  using 
Eq.  (3)  to  obtain 
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Now  define  the  collimator  correction  factor  for  the  energy  Ei   as 


DRfl(ya) 
R,  =  ^ .  (5) 

6  Th^^^ 


R   actually  represents  the  ratio  of  the  dosimeter  response  for  direct  radia- 
tion plus  that  radiation  reaching  the  detector  as  a  result  of  collimator 
interactions  to  the  response  resulting  from  direct  radiation  alone.   Com- 
bining Eq.  (4)  and  Eq.  (5)  yields 


Next  consider  a  point  source  which  emits  gamma  rays  with  an  energy 
distribution  identical  to  that  which  is  to  be  investigated  in  the  measure- 
ments for  determining  the  accuracy  of  the  unfolding  technique.   Assume  that 
this  point  source,  when  used  in  the  same  geometry  as  the  monoenergetic 
point  source  used  in  determining  Eq.  (3) ,  results  in  an  exposure  dose 

V^,   (mr)  such  that 
in 


P^^(mr)  =  P^^(mr)  .  (7) 


Calibrating  the  dosimeter  under  the  same  conditions  as  used  in  obtaining 
Eq.  (3)  results  in  the  relation 


PR^(vja)  =  k     P^,  (mr)  (8) 

m        sin 

where    V^    (ua)  =  dosimeter  response  in  microamperes  resulting  from  gamma 
m 

rays,  with  the  energy  distribution  previously  described, 
incident  normal  to  the  axis  of  the  dosimeter. 
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k     =   the  proportionality  constant  relating  the  dosimeter 
response  to  the  exposure  dose  for  this  point  source. 

Now  assume  that  a  calibrated  plane  source  is  available  which  emits  gamma 
rays  with  an  energy  distribution  identical  to  that  simulated  for  the  dosim- 
eter and  collimator-spectrometer  measurements.   Further  assume  that  this 
plane  source  is  such  that 

D^^(mr)  =  D^^(mr)  (9) 

E 1  s 

where  D^^(mr)  and  D^,  (mr)  are  the  theoretical  exposure  doses  which  would  be 
in  in 

expected  respectively  for  the  dosimeter  located  within  the  collimator  and 
exposed  to  the  monoenergetic  and  polyenergetic  plane  sources. 

g 

Define  DR  (ya)  as  the  response  of  the  dosimeter  located  within  the 

o  m 

collimator  resulting  from  an  exposure  to  the  polyenergetic  plane  source 
where  the  radiation  is  incident  at  an  angle  0  to  the  axis  of  the  dosimeter. 
Then  in  a  manner  similar  to  that  used  for  the  monoenergetic  plane  source 
the  relation  betx^/een  the  dosimeter  response  and  the  theoretical  exposure 

g 

dose,   D  ,  (pa),  can  be  represented  by 
W  in 

^DR^(,a)  =  R^  ^D^^<.a)  ,  (10) 


or 


k     D;  (mr) 

.DR^ya)  =  R^  -^~^ (11) 

3   m        s      s 

a„ 


where    R  =  the  collimator  correction  factor  for  the  gamma-ray  energy 
distribution  now  being  considered  (R   is  analogous  to  R 

S  II  1 


as  defined  by  Eq.  (5)), 
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OL     =   the  angular  correction  factor  of  the  dosimeter  for  gamma 

radiation  xi/ith  the  energy  distribution  now  being  considered, 

Using  Eq.  (9),  then  Eq.  (11)  can  be  written  as 


.DR:(ya) 

7   m 


R^  k     D;l(mr) 
s   s   Th 


(12) 


and  using  Eq.  (6)  a  substitution  for  D  /^(mr)  can  be  made  to  obtain 


R  k 

.DR^ya)  =  -^ 
3   m         8 


"-e'  e<'^^^^ 


El   El 


(13) 


Combining  Eqs.  (12)  and  (13)  yields 


,DR°(ya) 

i     m 


R 
s 

r  En 

s 

[\\ 

s 

,DR^l(ya) 

J   m 


R  fe    „ 
s   s  ^Ei ,   V 


a 


0 


(14) 


Hence,  using  Eq.  (14),  a  calibration  curve  could  be  obtained  which 
would  be  applicable  for  the  gamma-ray  energy  distribution  vjhich  is  to  be 
measured.   This  can  be  done  by  plotting  the  middle  group  of  terms,  i.e., 


,El 


El 


the  group  containing  DR  Hya) ,  versus  the  theoretical  dose,  D  Mmr) ,  and 


6  m 


Th 


then  using  a  regression  line  to  fit  these  points.   The  slope  of  this  line 
would  then  be  equivalent  to  (R  fe  )/ct„  . 

s   S     D 

Hov'jever ,  if  a  calibration  curve  is  to  be  determined  using  Eq.  (14), 

evaluation  of  the  ratios  (R  /R„  ),  (a.Vot.)  and  (fe  /fe_  )  is  necessary. 

sh^     yy       sbi 

El   s 
First  consider  (a^/a  ).   With  the  definition  of  the  two  terms  in  this  ratio, 

this  ratio  can  be  expressed  as 
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[ 5 ] 

-^  =  -^^ .  (15) 

a^      PR^(ya) 
e    [_^ J 


I'R^Cya) 


m 


Using  the  relations  for  k     and  fe   given  by  Eqs.  (3)  and  (8)  yields 

S         Hi 


[-T ] 


[—5 1 


Since  P„,  (mr)  =  C_/(mr),  Eq.  (16)  can  be  reduced  to 
in        in 


k  m^(ya) 

s      m 


(16) 


^Ei   PR^l(ya) 
m 


Thus  the  product  of  Eqs.  (15)  and  (17)  results  in 


(17) 


OS    D  m 


(18) 


a^  fe_    M^Hva) 
6   Ej    em 


Substitution  of  Eq.  (18)  into  Eq.  (14)  yields 


R  .PR^(ua)  R     k 

DR^(ya)  =  -^  [^-| ]  .Dr";!  (ya)  =  -^^-  D^^dnr)  .       (19) 

6   m       R     od^i  .   V  6   m  s     Th 

El  ei?R^^(ua)  Uq 


S         El 
In  order  to  evaluate  the  ratios  (R  /R^  )  and  („PR  (ya)/„'PR  Mya)),  it 

s   El       b  m     o  m 

is  necessary  to  know  the  energy  distribution  of  gamma  rays  emitted  from  the 
polyenergetic  plane  source.   It  was  noted  previously  that  it  would  be 
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possible,  although  quite  tedious,  to  obtain  this  energy  distribution  from 
theoretical  considerations.   However,  since  the  dosimeter  measures  a  quantity 
which  is  proportional  to  exposure  dose  integrated  over  time  and  energy,  only 
the  average  energy  of  the  energy  distribution  under  consideration  is  necessary, 
This  average  energy  can  be  approximated  by  the  single  scattered  gamma-ray 
energy  which  can  be  computed  using  the  scattering  angle  defined  by  the  experi- 
mental configuration  under  investigation.   The  validity  of  this  approximation 
will  be  discussed  in  Section  5.0.   Then  if  both  the  collimator  correction 
factor  and  the  dosimeter  response  can  be  expressed,  either  experimentally 

or  theoretically,  as  a  function  of  gamma-ray  energy,  the  ratios  (R  /R   )  and 

S   E^ 

s        El 
(-PR  (\ia)/J)R   •'■(ya))  can  be  evaluated  using  the  single  scattered  energies. 

y   m      a   m 

To  obtain  a  better  interpretation  of  the  physical  significance  of 
Eq.  (19),  consider  the  follov/ing  graphical  approach.   First  v;rite  Eq.  (6) 
in  the  form 


gDR^'(pa)  =k^^  D^'(mr)  (20) 

El 
where  k   =  R   k      fo.^    .      In  a  similar  manner,  Eq.  (12)  can  be  written  in 
tp]^     E]^   E  ]^   a 

the  form 


nR^(ya)  =  k  D^l(mr)  (21) 

9  m        s   in 

where  k  =  R  fe  /a„.   Solving  Eqs.  (20)  and  (21)  for  D^^(mr)  yields  the 
s    s   s   y  ih 

relation 


DR^l(ya)    .DR^ya) 

-- kl =  — k ^  ^Ti^^^^>  '  (22) 

El  s 

k 

or  .DR^(ya)  =  r--  .^^1  (ya)  .  (23) 

6   m       k^   D   m 
El 


Comparing  this  result  with 


17 


R    .PR"(ya) 
gDR^(ya)  ^  --   [--^— ^  e^m 


]  ,DR:l(ya) 


(24) 


from  Eq.  (19)  results  in 


R    .PR:(ya) 
k  =  ^  [^ ]  k^, 

D  m 


(2.5) 


Equation  (25)  expresses  the  relationship  betv/een  the  slopes  of  the  calibra- 
tion lines  for  the  monoenergetic  and  polyenergetic  plane  sources. 

Now  if  the  calibration  curves  represented  by  Eqs.  (20)  and  (21)  are 
plotted  on  the  same  graph  they  would  appear  as  illustrated  in  Figure  1.   It 
should  be  noted  that  k   is  the  slope  corresponding  to  the  dosimeter  calibra- 
tion  curve  for  the  monoenergetic  plane  source  and  k   is  the  slope  correspond- 
ing to  the  dosimeter  calibration  curve  for  the  polyenergetic  plane  source. 


1     m  a 


pe  =  k 


Figure  1.   Dosimeter  response  versus  exposure  dose  for  exposure  of  a 

dosimeter  to  a  monoenergetic  plane  source  (slope  =  k   )  and 
a  polyenergetic  plane  source  (slope  =  k  ).  ^ 
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In  representing  these  two  curves  graphically,  the  author  has  required  that 

s        El 
the  product  of  the  ratios  (R  /R^  )  and  (^Pr  (ya)/^PR  MPa))  be  less  than 

s   Li       D  m     o  m 

unity  as  this  is  the  case  which  v;ill  be  encountered  in  this  work. 

Examination  of  this  figure  indicates  that  the  two  calibration  curves 

are  of  the  same  form  and  in  fact  that  they  differ  only  with  respect  to  their 

slopes.   Referring  to  Eq.  (25)  it  is  apparent  that  the  only  difference  in 

the  slopes  is  that  to  obtain  k  from  k^   a  correction  is  necessary  to  account 
-  s       El 

for  the  dependence  of  the  collimator  correction  factor  and  the  dosimeter 
response  on  the  gamma-ray  energy. 

Therefore,  using  Eq.  (19)  and  a  monoenergetic  plane  source,  it  is 
possible  to  obtain  a  calibration  curve  for  the  dosimeter  which  will  be  suit- 
able for  use  in  the  experiment  described  in  Section  2.2.   This  calibration 
curve  v/ill  correct  for  the  energy  response  of  the  dosimeter  as  well  as  the 
effects  of  the  collimator  and  thus  will  make  possible  the  measurement  of  an 
exposure  dose  with  the  dosimeter  which  can  be  compared  to  that  computed 
using  the  collimator-spectrometer  system  measurements  and  the  unfolding 
technique. 
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3.0   EXPERIMENTAL  APPARATUS 

Experimental  measurements  were  conducted  at  the  Kansas  State  University 
Nuclear  Engineering  Shielding  Facility  on  a  3500  square  foot  concrete  slab. 
This  slab  was  equipped  with  two  fifty-foot  towers  located  approximately 
sixty  feet  apart  on  one  diagonal  of  the  concrete  slab. 

A  government  surplus  antiaircraft  searchlight  which  had  been  modified 
by  removing  the  searchlight  and  replacing  it  by  a  3500  pound  lead  collima- 
tor was  used  in  the  experimental  measurements.   The  collimator  housing  was 
of  such  a  design  as  to  enable  it  to  be  used  with  either  a  blank  collimator 
(no  aperture)  for  background  measurements  or  a  two  inch  diameter  by  fourteen 
and  one-half  inch  long  collimator  for  differential  dose  measurements.   This 
was  possible  through  the  insertion  of  the  appropriate  collimator  into  the 
collimator  housing.   The  instrument  cavity,  located  within  the  collimator, 
was  in  the  form  of  a  right  circular  cylinder  with  dimensions  of  three 
inches  by  thirteen  and  one-half  inches.   Copper  tubing  encircled  the 
instrument  cavity  to  allow  preheated  water  to  be  forced  through  the  colli- 
mator housing  walls  and  thereby  maintain  the  instrument  cavity  at  a  constant 
temperature.   A  Precision  Scientific  Company  constant  temperature  circu- 
lating system  was  used  to  heat  and  force  water  at  an  elevated  temperature 
through  the  copper  tubing.   A  working  drawing  of  the  collimator  system  is 
shown  in  Figure  2. 

A  Harshaw  integral  line  assembly,  Type  12S12,  was  used  in  measuring 
the  complex  gamma-ray  energy  spectra.   This  assembly,  consisting  of  a  three 
inch  by  three  inch  cylindrical  thallium-activated  sodium-iodide  crystal, 
photomultiplier  tube,  and  associated  parts,  was  used  with  a  Hamner  pre- 
amplifier. Model  N-356.   Power  was  supplied  to  the  photomultiplier  tube  by 
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a  John  Fluke  Model  412A  high  voltage  power  supply.   The  preamplifier, 
which  was  driven  by  a  Model  HB2AM  Kepco  voltage  regulated  power  supply, 
was  capable  of  sending  pulses  from  the  photomultiplier  tube  through  one 
thousand  feet  of  signal  cable.   Pulses  received  from  the  preamplifier  were 
analyzed  using  a  Technical  Measurements  Corporation  (TMC)  multichannel 
analyzer  system.   This  system  consisted  of  a  TMC  Model  402-6  400  channel  pulse- 
height  analyzer  and  a  TMC  Model  520  type-punch-read  control  unit.   The  data 
accumulated  by  the  multichannel  analyzer  system  were  read  out  using  both  a 
TMC  Model  530  computer  typewriter  and  a  TMC  540R  paper  tape  punch  in  con- 
junction with  the  type-punch-read  control  unit. 

Dosimeter  measurements  were  made  using  a  Victoreen  Instrument  Company 
air  equivalent  stray  radiation  chamber.  Model  239.   This  indirect  reading 
dosimeter  had  a  ten  milliroentgen  total  dose  capability  and  was  read  using 
a  Technical  Operations  Incorporated  dosimeter  charger-reader.   The  inte- 
grated dose  accumulated  by  the  dosimeter  was  read  on  the  charger-reader  in 
terms  of  microamperes. 

Two  cobalt-60  sources  of  different  strengths  were  used  in  making  the 
required  measurements.   The  use  of  two  sources  was  necessary  because  of  the 
differences  in  sensitivity  of  the  dosimeter  and  the  scintillation  detector. 

Measurements  made  with  the  scintillation  detector  used  a  Technical 
Operations  Incorporated  gamma-ray  projector,  Model  402,  with  a  nominally 
0.3  curie  cobalt-60  source.   (See  reference  11.)   This  unit  was  operated  by 
cranking  the  source  out  of  its  lead  storage  container  by  means  of  a  flexible 
cable  on  the  end  of  which  the  source  was  mounted.   The  flexible  metallic 
hose  through  which  the  source  was  moved  was  equipped  with  an  electrical 
control  system  (battery  powered)  so  that  the  operator  cculd  determine  when 
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the  source  was  out  of  the  storage  container  and  when  it  was  in  position  at 
the  end  of  the  metallic  hose. 

Dosimeter  measurements  were  made  using  a  Technical  Operations  Incor- 
porated duplex  hydraulic  source  circulation  system,  Model  539.   In  this 
system  water  was  drawn  from  a  fifty-five  gallon  reservoir  through  a  high 
speed  pump  thereby  applying  pressure  to  the  rear  of  a  source  piston  assembly. 
A  nominally  ninety  curie  cobalt-60  source  was  mounted  on  this  assembly. 
Thus  the  source  was  forced  from  its  storage  container  into  a  length  of 
polyethylene  tubing  which  had  been  positioned  to  satisfy  the  requirements 
of  the  experimental  configuration.   By  means  of  a  forward  and  reverse 
switch  the  source  could  be  positioned  to  within  one  inch  of  a  selected 
location  in  the  tubing.   A  detailed  description  of  this  source  circulation 
system,  with  operating  instructions,  is  contained  in  reference  12. 

The  calibration  of  the  dosimeter  in  the  collimator  used  a  simulated 
plane  source  of  cobalt-60.   A  uniform,  plane  isotropic  source  of  cobalt-60 
was  approximated  by  traversing  a  "line"  source  over  a  rectangular  plane. 
The  "line"  source  consisted  of  eighteen  nominally  four  millicurie  cobalt-60 
sources.   These  sources  were  positioned  on  a  cardboard  holder  by  placing 
the  small  wire  loop  attached  to  the  end  of  each  source  over  one  of  the 
eighteen  wire  hooks  located  at  one-half  inch  intervals  along  the  length 
of  the  holder.   This  line  source,  which  is  illustrated  in  Figure  3,  was 
then  attached  to  the  plane  simulator  shown  in  Figure  4  which  traversed  the 
line  source  in  the  horizontal  and  vertical  directions  simultaneously.   The 
vertical  speed  was  approximately  one  inch  per  minute  and  one  complete 
revolution  of  the  "line"  source  in  the  horizontal  direction  took  about 
nineteen  seconds.   The  plane  simulator  was  capable  of  simulating  an  approxi- 
mately twenty-eight  inch  square  plane  isotropic  source  of  cobalt-60. 
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Figure  3.   'Line'  source  consisting  of  eighteen  cobalt-60  sources 


With  regard  to  radiation  safety  ,  each  person  in  the  experimental  area 
was  equipped  with  a  film  badge,  two  pocket  dosimeters  and  a  radiation 
survey  m.eter.   Guards  were  employed  to  prevent  the  entry  of  unauthorized 
personnel  to  the  experimental  area.   This  was  accomplished  by  having  the 
guards  constantly  monitor  the  tv/o  milliroentgen  per  hour  perimeter  surround- 
ing the  experimental  area  while  the  experiments  were  in  progress.   In 
addition  all  experiments  were  carried  out  in  accordance  with  the  rules  and 
regulations  set  forth  by  the  Kansas  State  University  Radiation  Safety 
Committee  and  the  appropriate  state  and  federal  licenses. 
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4.0   EXPERIMENTAL  PROCEDURE 

Polyenergetic  gamma  rays  were  obtained  by  scattering  cobalt-60  radia- 
tion from  a  concrete  slab.   Two  fifty  foot  towers  were  used  to  suspend  a 
cobalt-60  source  above  the  concrete  slab.   For  selected  scattering  angles 
(see  Figure  5)  the  collimator   defined   corresponding  areas  on  the  slab. 

Experimental  measurements  were  made  using  both  the  scintillation 
detector  and  the  ten-milliroentgen-f ull-scale  dosimeter.   These  detection 
instruments  were  positioned,  one  at  a  time,  in  the  instrument  cavity  of  the 
collimator  and  measurements  were  made  for  each  of  the  experimental  con- 
figurations to  be  considered.   The  geometry  of  the  experimental  configura- 
tions is  illustrated  in  Figure  5.   Each  of  these  configurations  was  defined 
by  the  angles  6q,  6  and  (^  ^.s   defined  in  this  figure. 

It  should  be  noted  that  because  of  the  design  of  the  collimator  the 
detector  height  for  all  measurements  was  maintained  at  five  feet  two  inches, 
Also  the  source  height  above  the  concrete  slab  was  held  constant  at  ten 
feet. 

For  each  experimental  configuration  and  for  both  the  scintillation 
detector  and  the  dosimeter,  a  minimum  of  six  measurements,  three  Xi7ith  the 
two  inch  collimator,  and  three  using  the  blank  collimator,  were  made  to 
obtain  a  measure  of  the  dose  and  background  for  the  configuration  under 
investigation.   The  geom.etry  of  each  experimental  configuration  was  repro- 
duced independently  at  least  three  times  for  both  the  scintillation 
detector  and  dosimeter  measurements.   This  was  accomplished  by  making 
only  one  dose  and  background  measurement  for  a  particular  experimental 
configuration  before  proceeding  to  the  next  one. 
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Figure  5.   Experimental  configuration  used  in  the  unfolding  confirmation, 
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Spectral  measurements  were  made  using  the  sodium  iodide  scintillation 
detector  and  the  Technical  Operations  gamma-ray  projector  with  the  nominally 
three-tenths  curie  cobalt-60  source.   The  source  tip  of  the  metallic  hose 
of  this  unit  was  held  in  place  above  the  concrete  slab  through  the  use  of  a 
rope  from  each  of  the  fifty  foot  towers  and  guide  wires  to  the  concrete 
slab.   For  each  measurement  the  collimator  was  positioned  and  the  source 
was  cranked  from  its  storage  container  through  the  metallic  hose  to  its 
position.   Following  the  completion  of  the  experimental  run  the  source  was 
returned  to  its  storage  container  and  preparations  were  made  for  the  next 
measurement.   A  quality  control  program  was  developed  to  assure  a  consistent 
level  of  operation  for  the  gamma-ray  scintillation  spectrometer  system. 
This  quality  control  program  is  presented  in  Appendix  A. 

The  ten-milliroentgen  full-scale  dosimeter  was  calibrated  using  the 
"line"  source  and  the  plane  simulator.   The  dosimeter  was  first  positioned 
in  the  instrument  cavity  of  the  collimator  so  that  the  center  of  its  active 
volume  coincided  with  the  axis  of  rotation  of  the  collimator  and  its  verti- 
cal axis  coincided  with  the  collimator  axis.   The  calibration  then  consisted 
of  measuring  the  dosimeter  response  for  an  exposure  of  one,  two,  three  and 
four  traversals  of  the  line  source  over  the  prescribed  plane.   Figure  6 
illustrates  the  experimental  configuration  for  dosimeter  calibration. 

Measurements  made  with  the  ten-milliroentgen  full-scale  dosimeters 
were  carried  out  using  the  Technical  Operations  Incorporated  duplex  hydraulic 
source  circulation  system  with  the  nominally  ninety  curie  cobalt-60  source. 
For  these  measurements  the  dosim.eter  was  charged  and  put  inside  the  colli- 
mator and  then  the  collimator  was  positioned.   The  source  was  then  pumped 
from  its  storage  container  into  the  polyethylene  tubing  and  positioned 
using  the  forward  and  reverse  switch  on  the  pump  console  and  a  spotting 
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Figure  6.   Experimental  configuration 

for  the  dosimeter  calibration 


telescope.   The  tubing  was  positioned  above  the  concrete  slab  by  attaching 
it  to  a  heavy  vjire  support  which  was  in  turn  suspended  in  air  through  the 
use  of  light  ropes  from  the  two  fifty  foot  towers.   Guide  wires  from  the 
wire  support  to  the  concrete  slab  were  used  so  that  the  tubing  could  be 
kept  in  a  fixed  position.   At  the  end  of  a  given  experimental  run  the  source 
was  returned  to  its  storage  container  and  the  dosimeter  was  then  removed 
from  the  collimator  and  read  using  the  Technical  Operations  Incorporated 
charger-reader . 

Before  the  dosimeter  calibration  equations  could  be  determined 
it  was  first  necessary  to  calibrate  the  cobalt-60  simulated  plane 
source  in  the  manner  described  in  Appendix  B.   This  determination  of  the 
cobalt-60  simulated  plane  source  strength  required  that  first  the  dosimeter 
be  calibrated  using  a  calibrated  cobalt-60  point  source  (Appendix  C) . 
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Since  this  work  contains  a  number  of  different  experiments,  some  of  which 
are  dependent  on  others,  the  experiments  and  calculations  performed  during 
the  course  of  this  work  are  illustrated  in  schematic  form  in  Figure  7. 
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5.0   CALIBRATION  OF  THE  DOSIMETER 
LOCATED  WITHIN  THE  COLLIMATOR 

The  ten-milliroentgen  total-dose  dosimeter  was  calibrated  with  the 
dosimeter  located  at  the  center  of  the  collimator.   As  a  result  of  the 
calibration  measurements  made,  equations  for  the  calibration  curves  could 
be  computed  which  would  be  suitable  for  use  in  the  dosimeter  measurements 
of  the  scattered  gamma-ray  energy  distributions  considered  in  Section  1.2. 
These  calibration  curves  are  based  on  Eq.  (19),  i.e., 


R 
,DR^(ya)  =  ^ 


m'"  '    R„ 

^1 


V      m 


,^R^^(ya) 
3      m 


.DR^jlCya)  =  -~-^  D!;l(mr) 
o  m  s    ih 


which  was  derived  in  Section  1.3. 

This  calibration  experiment  consisted  of  measuring  the  response  of  the 
dosimeter,  located  within  the  collimator,  exposed  to  the  simulated  plane 
isotropic  source  of  cobalt-60.   The  geometry  for  these  measurements  is 
illustrated  in  Figure  6  and  the  "line"  source  and  plane  simulator  used  in 
simulating  the  plane  isotropic  source  of  cobalt-6G  are  illustrated  in 
Figures  3  and  A  respectively.   Dosimeter  responses  were  measured  for  one, 
two,  three  and  four  simulations  of  the  plane  source.   These  measured  dosim- 
eter responses,  in  microamperes,  were  corrected  for  temperature  and 
pressure  using  Eq.  (C-1)  of  Appendix  C  and  normalized  to  a  twenty-eight 

minute  simulation  of  the  plane  source.   These  corrected  dosimeter  responses 

p 

correspond  to  .DR  ^ya)  in  Eq.  (19). 

D  m 

El 
The  theoretical  dose,  D  /^(mr)  in  Eq.  (19).  corresponding  to  the  experi- 
mental measurements  described  above  was  based  on  the  dose  which  would  be 
expected  at  the  center  of  the  dosimeter  due  to  one  simulation  of  the  plane 


32 


source  of  cobalt-60.   Thus  using  the  strength  of  the  simulated  plane  source 
as  determined  in  Appendix  B  the  theoretical  dose  can  be  calculated  from 


D^ii(mr)  . 


D 

d^A  ,  (26) 


Airr 


where    D  =  strength  of  the  plane  source  in  milliroentgens  per  twenty- 
eight  minutes, 
A  =  area,  in  square  centimeters,  on  the  plane  source  defined  by 

the  solid  angle  of  the  collimator, 
r  =  the  distance  in  centimeters  from  the  differential  area  d^A 
to  the  center  of  the  collimator. 
Recognizing  that  d^A  =  2'JTpdp,  r  =  z  sec  6  and  p  =  z  tan  6,  then  a  change 
of  variable  from  p  to  6  in  Eq.  (26)  yields 


D^^(mr)  = 


1/2  D   tane  dG  =  1/2  D  (-ln(cos9,))  (27) 

P  P         ^ 

0 


in  which  6, ,  p  and  z  are  defined  in  Figure  6.   Substitution  of  the  appro- 
priate value  for  0^  yields 


D^^(mr)  =  8.66  x  10""^  D  (28) 

in  p 


or 


Df;^(mr)  =  0.568  +  0.014  (29) 

ih  — 


Although  Dp,  v^7hich  could  not  be  measured  experimentally,  has  the  units  of 
dose,  it  actually  represents  the  product  of  source  strength,  exposure  time 
and  the  flux  to  dose  conversion  factor.   This  unusual  term  was  employed 
for  convenience  in  carrying  out  the  calculations. 
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where  the  value  for  D  was  obtained  from  the  calibration  of  the  plane  source 

P 

(see  Appendix  B) .   Thus  Eq.  (29)  represents  the  expected  dose  per  simulation 

of  the  plane  source  due  only  to  those  gamma-rays  x-;hich  leave  the  plane 

source  in  such  a  direction  as  to  strike  the  detector  without  the  need  of 

an  intermediate  scatter. 

It  was  noted  in  Section  2.3  that  the  evaluation  of  the  ratios  (R  /R^  ) 

s  El 

s        El 
and  (q1^  (pa)/  PR  ^(pa))  requires  a  knowledge  of  the  scattered  gamma-ray 
y  m     0  m 

energy  distributions  for  which  the  calibration  curves  are  to  be  applicable. 
Also  in  this  previous  discussion,  it  was  stated  that  since  only  the  average 
energy  of  the  gamma-ray  energy  distribution  was  necessary  for  the  evaluation 
of  these  two  ratios,  the  single  scattered  energy  based  on  the  experimental 
geometry  could  be  used  to  approximate  the  average  energy  of  each  experi- 
mental configuration.   However,  in  evaluating  this  approximation  it  is  first 
necessary  to  consider  how  the  dosimeter  response  and  the  collimator  correc- 
tion factor  vary  with  gamma-ray  energy. 

The  problem  of  collimator  effects  has  been  examined  on  a  theoretical 
basis  by  Mather  (13)  and  experimentally  by  Tomnovec  and  Mather  (14).   In 
both  of  these  papers  results  have  been  given  relating  the  ratio  of  collimator- 
interacted  radiation  to  direct  radiation  (this  corresponds  to  R  -  1.0  as 
defined  in  this  report)  at  some  point  within  the  collimator  to  the  energy 
of  the  gamma  rays  incident  on  the  face  of  the  collimator.   From  these 
studies  it  can  be  seen  that  the  collimator  correction  factor  R  can  be  very 
closely  approximated  as  a  linear  function  of  gamma-ray  energy.   This  approx- 
imation will  not  introduce  an  error  of  more  than  five  percent  in  the  value 
.s,   ,  ,  ^„Ei 


of  (  PR^ya)/J?R';l(ya))  ->  R  /R^ 
D   m     9   m  s   Ej 
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If  R  is  assumed  to  be  a  linear  function  of  gamma-ray  energy,  all  that 
is  needed  to  relate  R  and  the  gamma-ray  energy  are  two  points.   One  point 
on  a  curve  of  R  versus  gamma-ray  energy  is  the  point  (1.000,  0.00).   This 
point  follows  directly  from  the  definition  of  R  and  from  the  fact  that  the 
absorption  cross  section  for  lead  at  low  energies  greatly  outweighs  the 
Compton  cross  section  and  that  the  linear  attenuation  coefficient  increases 

rapidly  for  decreasing  gamma-ray  energy.   The  value  of  R   can  be  calculated 

El 

from  the  experimental  measurements  made  using  the  plane  source  and  knowing 

the  strength  of  the  simulated  plane  source.   Hence  it  is  found  that 

R   =  1.677  and  therefore  the  second  point  is  (1.677,  1.25)  where  1.25  is 
El 

the  average  energy  in  MeV  of  the  two  gamma  rays  emitted  by  cobalt-60. 
Using  these  two  points  then  the  slope  of  the  line  can  be  computed  and  the 
resulting  expression  relating  the  collimator  correction  factor  R  to  the 
incident  gamma-ray  energy  E  is 

R  =  1.000  +  0.542  E  .  (30) 

s        El 
The  ratio  iJ)R    (\ia)  /  J)R   ^    (ua))  can  be  obtained  directly  from  a  curve 
o   m     6   m 

(15)  of  dosimeter  response  at  an  energy  E,  relative  to  the  response  to  cobalt- 
60  radiation  versus  gamma-ray  energy.   This  curve  indicates  that  the  variation 
of  the  dosimeter  response  relative  to  its  response  to  cobalt-60  radiation  is 
nearly  constant  for  decreasing  gamma-ray  energies  down  to  approximately 
0.2  MeV. 

Based  on  this  analysis  of  the  energy  dependence  of  the  ratios  (R  /R   ) 

s        El 
and  (qPR  (ya)/.PR  ^    (pa)),  it  is  possible  to  predict  the  error  introduced 

u   m     o  m 

through  the  use  of  the  single  scattered  energy  as  the  average  energy  of  the 
gamma-ray  energy  distribution  under  investigation.   Since  the  ratio 
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S  cJl 

(  PR  (ija)/PR  Hpa))  is  essentially  a  constant  for  gamma-ray  energies  down 

cm     o   m 

to  approximately  0.2  MeV,  then  as  long  as  the  average  energy  is  sufficiently 
larger  than  this  lower  limit  an  error  incurred  through  the  use  of  the  single 
scattered  energy  will  affect  this  ratio  negligibly.   The  ratio  (R  /R   ), 

S    i-i  1 

however,  is  sensitive  to  an  error  in  the  average  energy  used.  As  an  example 
if  the  single  scattered  energy  calculated  is  a  factor  of  two  larger  than  the 
true  average  energy  of  the  scattered  gamma-ray  energy  distribution  the  error 
introduced  in  the  ratio  (R  /R„  )  will  be  approximately  ten  percent. 

S    ■'-'1 

The  single  scattered  energy  approximation  is  an  upper  limit  to  the 
average  energy  of  the  gamma  rays  reaching  the  detector  after  scattering  from 
the  concrete.   This  follows  from  the  fact  that  a  significant  fraction  of  the 
gamma  rays  striking  the  detector  have  suffered  multiple  collisions  in  the 
concrete.   These  multiply  scattered  gamma  rays  will  have  a  much  lower  energy 
and  thus  it  is  obvious  that  the  true  average  gamma-ray  energy  of  the  concrete 
scattered  radiation  will  be  less  than  the  single  scattered  energy  and  that 
the  single  scattered  energy  will  be  an  upper  limit. 

Nov?  the  single  scattered  energy  E   ,  for  each  of  the  experimental 
configurations  being  considered,  can  be  computed  from 


^1 
^ss  "       Eid.O  -  cos  e  )  ^^-^^ 

1.0  +  2 — 

mo  c^ 

where     E^  =  incident  gamma-ray  energy  on  the  concrete,  1.25  MeV  average 
for  cobalt-60, 
mQ  c^   =  rest  energy  of  an  electron,  0.51  MeV,  and 
cos  9   =  cosine  of  the  scattering  angle  6  . 
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The  cosine  of  the  scattering  angle  can  be  computed  using 

cos  e   =  -  sin  Bq  sin  6  cos  (^  +  cos   9^  cos  9  (32) 

where  the  angles  Gq,  9  and  <^   are  defined  in  Figure  5. 

Thus  the  single  scattered  gamma-ray  energy  for  each  of  the  experimental 
configurations  can  be  determined  using  Eqs.  (31)  and  (32).   The  experimental 
configurations,  designated  as  P^,  P2  and  P3,  and  their  defining  angles  9  , 
9  and  (j)  along  with  the  single  scattered  energies  for  each  are  listed  in 
Table  I. 


Table  I.   Experimental  configurations  and  corresponding 
energies  of  once  scattered  gamma  rays 


.    ^  ,  Defining  angles  in  degrees  E 

Experimental  000  gg 

Configuration  9q  9  (^  (MeV) 


?!  00.0       72.0       00.0  0.295 

?2  ^1-4       72.0       00.0  0.505 

P3  60.0       72.0       52.5  0.481 


Using  the  single  scattered  energies  as  presented  in  Table  I,  the  values 

s        El 
of  (.PR  (ya)/-PR  ^(ya))  and  R  can  be  determined  in  the  manner  previously 
6   m     9   m  s 

described.   Hence  with  the  value  of  R  and  the  previously  calculated  value 

S  r  J 

of  R^   the  ratio  (R  /R^  )  can  be  evaluated  for  each  of  the  experimental 
E^  s   E] 

configurations.   The  values  of  these  ratios  for  each  of  the  experimental 
configurations  are  presented  in  Table  II. 
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Table  II.   Values  of  (R  /R^  )  and  (^PR^(ya) /^PR^^ (ya)) 

s   Ej       D   m     e   m 

for  each  experimental  configuration  based 

on  energies  of  once  scattered  gamma  rays 


R                     ^PR^(ya) 
Experimental  s  6  m 


configuration  R„  -i-j-Ei,   . 

El  PR  Mya) 

u  m 


Pi  0.692  0.993 

P2  0.759  0.989 

P3  0.752  0.992 


With  the  values  given  in  the  above  table,  the  experimentally  measured 

El  El 

dosimeter  responses,  „DR  ^  (ya) ,  and  the  theoretical  doses,  D  Mmr)  ,  from 

o  m  in 

Eq.  (29),  the  equation  of  the  calibration  curve  can  be  calculated  for  each 
of  the  three  experimental  configurations.   This  is  accomplished  through  the 

application  of  the  regression  analysis  of  Appendix  D  to  Eq.  (19).   This 

El 
application  involves  considering  x  as  D  /-(mr)  and  y  as  the  product  (R  /R   ) 

X  n  s   Hi  1 

(  pR^ya)/  pR;^l(ya))  .DR^^ya). 
D  m     y  m       o  m 

As  an  example,  the  determination  of  the  equation  for  the  calibration 

curve  for  the  P^  experimental  configuration  is  presented  below.   The 

El 
experimental  dosimeter  responses  .DR  ^  (ya) ,  corrected  for  temperature  and 

y  m 

pressure,  the  product  (R  /R^  )  (^PR^  (ya) /^PR^^  (ya) )  „DR^^  (ya)  and  Di?i(mr) 

sEiDmym       ym  in 

are  listed  in  Table  III  for  one,  two,  three  and  four  simulations  of  the  plane 
source  of  cobalt-60.   The  last  two  columns  in  this  table  represent  y  and  x 
respectively  as  used  in  the  regression  analysis  which  follows. 
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Table  III.   Data  used  to  determine  the  equation  of  the 
calibration  line  for  the  Pi  configuration 


g 

Simulations    of          T^T.El  /      \            .    s    ..6      m             .      T^T^E^  ,      .  txEi  ,      . 

1,1           c                      o^R      (^^)           (^ — )( i )  D^R      (pa)  Drru(mr) 

Plane   Source          6     m                       R^         ^-i^Ei,      v  6     m  Th 

El      ^\'(ya) 


8.40  5.78  0.568±  0.014 

7.83  5.38 

8.38  5.76 

8.23  5.66 

15.58  10.71  1.137±  0.028 

15.66  10.77 

15.55  10.69 

16.11  11.08 

22.63  15.56  1.705±  0.042 

23.49  16.15 

23.38  16.08 

23.70  16.30 

31.37  21.56  2.272±0.056 

31.54  21.69 

30.90  21.25 


Since  the  x  values  (D  Mmr))  in  Table  III  have  a  constant  percentage 
error,  the  equation  of  the  calibration  line  can  be  determined  using  the 
regression  analysis  for  a  controlled  experiment  as  outlined  in  Appendix  D. 
The  slope  of  the  line  is 

15 

y  X.  y. 

b  =  -^y =  9.46  .  (33) 

I      (Xi)2 
1=1 

The  estimate  of  o^,  from  Eq.  (D-20)  of  Appendix  D,  is 

(s')2  =  s2  +  b2  V[x]  (34) 
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or 

15  15  15 

I      (7^)2  -  [(  I   X.  y.)2  /  I      (x.)2] 

(s')2  =  ^ ^^ ^^ +  b2  V[x]       (35) 

£  -  1 


where  V[x]  is  the  variance  on  x  and  (£-1)  represents  the  degrees  of  freedom. 
Substitution  of  the  appropriate  values  yields 

s'  =  [0.0688  +  b2(0.0251  x/r*"^^  .  (36) 

However,  the  calibration  line  is  to  be  used  in  reverse  to  convert  a 
measurement  in  microamperes  to  an  exposure  dose  in  milliroentgens.   Therefore 
using  Eq.  (D-21)  of  Appendix  D,  the  value  of  x  can  be  determined  for  y',  the 
mean  of  m  new  observations.   The  form  of  this  equation  is 

""   b  -  b    U      15      ^  ^"^^^ 


b2  I       (x  ) 
i=l 


2 


if  g',  as  defined  by  Eq.  (D-22)  of  Appendix  D,  is  less  than  0.1.   Substitu- 
tion of  the  appropriate  values  into  the  expression  for  g'  indicate  that,  in 
the  range  of  interest,  g'  is  less  than  0.1.   Thus  Eq.  (37)  is  a  valid 
representation  for  the  regression  line  in  reverse.   Recognition  that  x  in 
Eq.  (36)  can  be  replaced  by  y'/b,  substitution  of  the  proper  values  into 
Eq.  (37)  yields 


^  -  ^lL  +  (1.03A)[0.0688+(0.0251  y')^]^^^  .1  ,    (y')^      /^^   .... 
""   9.46'-  9.46  ^m   (9.46)^(33.60)"^     ^^^^ 


where  the  value  of  t  for  68  percent  confidence  and  fourteen  degrees  of 
freedom  is  1.034. 
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In  the  interest  of  clarity  the  equation  of  the  calibration  line  can  be 

El 
written  in  terms  of  the  experimental  parameters  by  replacing  x  by  D  Mmr) 

and  y'  by  DR  (ya)  which  from  Eq.  (19)  is  equal  to  the  product  term 

(R  /R^  )  [VR^(ua)/r)R^^(]ja)]J)R^Hva).      Thus  Eq.  (38)  can  be  written  as 

sii]^fc)m  om  cm 


D!^l(mr)    =^-^ 

^^  9.46 


(1.034)  [0.0688+(0.0251   ^DR^(ya))2]^/^     1            (  nR^(ya))2      1/2 
+  _L_^ [_  +  L_iB J  (39) 

9.46  m        (9.46)2(33.60) 


where  DR  (ya)  is  the  mean  of  m  new  observations.   Equation  (39)  represents 

D   m 

the  calibration  line  for  the  P^  experimental  configuration.   By  replacing 

the  values  of(R  /R^  )  („PR^ (ya) /„PR^^ (ya))  „DR^^ (ya)  in  Table  III  by  their 
sEj^umom       6m 

corresponding  values  for  the  P2  and  P3  configurations,  the  equations  repre- 
senting these  calibration  lines  for  these  configurations  can  be  obtained  in 
a  similar  manner.   Hence  the  equation  of  the  calibration  line  for  the  ?2 
experimental  configuration  is 


T>Ei  /      \        0  m 

D  Mmr)  =  

10.33 


(1.034)  [0.0823+(0.0251   ^DR^  (ya))^]  ^^^      1               (   DR^(ya))2      1/2 
+  1_JE [_  +  6_iB ]  (40) 

10.33  m        (10.33)2(33.60) 

and  for  the  P3  configuration  the  equation  is 


D!^l(mr)  =^-^ 

^^         10.26 


(1.034) [0.0811+(0.0251  „DR^(ya))2] ^^^      1  (^DR^(ya))2      1/2 

D      m r ,   o      m 

10.26  m        (10.26)2(33.60) 


+ 6_^ [_  +  6-^5 ]         .       (41) 
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6.0  DISCUSSION  OF  RESULTS 

The  experimental  measurements  were  carried  out  using  the  experimental 
configuration   illustrated  in  Figure  5  and  the  procedure  outlined  in  Section 
4.0.   Three  different  configurations,  obtained  by  varying  9   and  (j)  were 
investigated.   The  three  experimental  configurations  examined  were  designated 
as  P^,  ?2   and  P3  and  correspond  to  combinations  of  the  angles  6  ,   6   and  <^ 
as  shown  in  Table  I.   Measurements  were  made  using  the  scintillation  detector 
system  and  the  ten  milliroentgen  dosimeter  alternately  in  the  lead  collimator. 
These  measurements  utilized  both  the  two-inch  and  the  blank  collimators. 

The  experimentally  observed  gamma-ray  energy  spectra  were  obtained  using 
the  scintillation  detector  and  the  nominally  ninety  curie  cobalt-60  source. 
A  consistent  level  of  operation  of  the  scintillation  detector  and  the  pulse- 
height  analyzer  was  obtained  through  the  use  of  the  quality  control  program 
presented  in  Appendix  A. 

Gamma-ray  spectra  obtained  for  each  of  the  experimental  configurations 
using  the  two-inch  collimator  were  averaged  and  normalized  to  an  exposure 
time  of  one  hour.   The  corresponding  blank  collimator  measurements  were  also 
averaged  and  normalized  to  a  one  hour  exposure  time.   The  two  averaged  spectra 
for  each  experimental  configuration  were  then  differenced  and  put  into 
twenty-five  energy  bins  of  forty-eight  keV  per  bin.   These  binned  spectra 
were  then  used  as  input  to  the  unfolding  codes  discussed  in  Section  1.1.   The 
output  from  the  unfolding  codes  was  the  unfolded  pulse-height  distribution 
N.   The  averaged  and  differenced  experimental  spectra  and  their  corresponding 
unfolded  pulse-height  distributions  are  illustrated  in  Figures  8,  9  and  10. 

With  reference  to  Figures  8,  9  and  10,  in  both  the  original  spectra  and 
the  unfolded  pulse-height  distributions,  three  peaks  are  clearly  discernable. 
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The  lowest  energy  peak  is  due  to  gamma  rays  which  were  multiple  scattered 
in  the  concrete  and  the  middle  peak  corresponds  to  the  single  scattered 
gamma  rays.   The  small  peak  at  an  energy  of  approximately  nine-tenths  of  an 
MeV  corresponds  to  those  gamma  rays  which  were  singly  scattered  from  the  lip 
of  the  collimator.   Upon  examination  it  is  found  that  there  is  close  agreement 
between  the  concrete  single  scattered  peak  energies  (energy  corresponding 
to  the  middle  peak  in  each  spectrum)  and  the  single  scattered  energies  pre- 
sented in  Table  I. 

It  should  be  noted  that  the  contribution  to  the  pulse-height  distribu- 
tion from  the  first  energy  bin  in  Figures  8,  9  and  10  is  negative.   Although 
these  negative  contributions,  which  correspond  to  negative  fluxes,  are  not 
consistent  with  physical  reality  they  must  be  considered  when  the  exposure 
doses  for  each  of  the  unfolded  pulse-height  distributions  are  computed. 
This  results  from  the  fact  that  the  objective  of  this  report  is  to  determine 
how  accurate  the  unfolding  technique  of  Baran  (10)  is  when  used  to  compute 
exposure  doses  rather  than  a  detailed  analysis  of  the  unfolding  technique. 

Even  though  a  detailed  analysis  of  the  unfolding  technique  is  not 
within  the  immediate  objectives  of  this  work,  an  indication  as  to  the  source 
of  this  error  is  in  order.   The  lowest  energy  source  used  in  the  calibration 
of  the  spectrometer  system  was  cobalt-60  metastable  which  emits  gamma  rays 
with  energies  of  58.9  keV.   Since  the  first  energy  bin  covers  the  energy 
range  from  0.0  to  48.0  keV,  then  the  development  of  the  response  matrix 
required  an  extrapolation  of  the  experimental  data  to  obtain  a  response  func- 
tion corresponding  to  this  first  energy  bin.   A  linear  extrapolation  was  used 
in  developing  the  response  matrix  for  the  unfolding  technique  considered  in 
this  work  (10).   However,  such  a  linear  extrapolation  does  not  appear  to  be 
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consistent  with  the  interaction  processes  of  low  energy  gamma  rays  with  the 
scintillation  detector  and  in  fact  such  a  linear  extrapolation  results  in  an 
overestimation  of  the  gamma-ray  contribution  in  the  response  function  corre- 
sponding to  the  first  energy  bin.   This  problem  is  being  considered  in  detail 
by  Baran  (16). 

The  unfolded  pulse  height  distributions  having  units  of  gamma  rays  per 
square  centimeter  per  hour  were  converted  to  exposure  doses  using 

25 

D  =  G   I  N.  E.  (y,(E,)/p)^,^  (42) 

1=1 

where     D  =  exposure  dose  in  mr/hr, 

G  =  conversion  factor,  1.91  x  10~   (mr-gm/MeV) , 

N.  =  number  of  gamma  rays  per  square  cnetimeter  per  hour  in 

the  ith  energy  bin  and  . 

^a^^i^ 
( )    =  gamma-ray  energy  absorption  mass  attenuation  coefficient 

air 

(cm^/gm)  corresponding  to  the  energy  E.. 

The  error  associated  with  each  of  these  exposure  doses  was  determined 
using  an  error  analysis  derived  by  Baran  (10).   This  analysis  considers 
both  the  errors  in  the  experimental  measurements  and  those  introduced  by  the 
experimentally  determined  response  matrix  as  well  as  the  error  introduced 
by  the  unfolding  process  itself. 

Dosimeter  measurements  were  made  by  replacing  the  scintillation  detector 
by  the  ten-milliroentgen  total-dose  dosimeter.   The  same  experimental  con- 
figurations were  investigated  in  these  measurements  as  were  examined  using 
the  scintillation  detector.   However,  because  of  the  relative  sensitivities 
of  the  two  detection  instrum.ents ,  the  nominally  ninety  curie  cobalt-60 
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source  was  employed  in  these  measurements.   A  minimum  of  three  measurements 
were  made  for  each  configuration  using  the  two-inch  collimator.   Following 
a  series  of  measurements  utilizing  the  blank  collimator  it  was  concluded  that 
the  background  was  negligible.   This  conclusion  was  further  verified  by 
theoretical  calculations. 

Experimentally  measured  dosimeter  responses,  in  microamperes,  were 
corrected  for  temperature  and  pressure  using  Eq.  (C-1)  of  Appendix  C  and 
then  averaged  for  each  of  the  experimental  configurations.   The  average 
response  was  then  substituted  into  the  appropriate  equation  of  the  calibra- 
tion line  (Eq.  39,  40  or  41)  corresponding  to  the  experimental  configuration 
under  consideration  to  obtain  the  exposure  dose.   These  doses,  one  for  each 
configuration,  were  then  normalized  to  a  one  hour  exposure  time. 

In  order  that  a  meaningful  comparison  could  be  made  between  the  exposure 
doses  measured  using  these  two  methods,  a  correction  was  necessary  to  account 
for  the  difference  in  the  strength  of  the  sources  used  in  the  scintillation 
detector  and  dosimeter  measurements.   This  adjustment  was  made  by  normalizing 
the  exposure  doses  to  a  source  strength  of  one  curie.   Thus  the  exposure 
doses  obtained  from  the  unfolding  codes  and  the  dosimeter  measurements  were 
divided  by  0.243  +  0.007  and  82.6  +  1.8  curies  respectively.   These  source 
strengths  were  obtained  by  correcting  the  calibrated  strengths  of  these 
sources  for  the  decay  in  their  strengths  which  had  occurred  from  the  date 
they  were  calibrated  to  the  date  when  they  were  used  in  making  these  measure- 
ments.  These  normalized  exposure  doses  in  mr  per  hour  per  curie,  are  pre- 
sented in  Table  IV  for  each  of  the  experimental  configurations  considered. 
In  this  table  the  doses  obtained  as  a  result  of  the  application  of  the 
unfolding  technique  to  the  complex  spectra  observed  are  designated  as  the 
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Table  IV.   Comparison  of  spectral  and  dosimeter  exposure  doses 


Experimental        Spectral  Dose       Dosimeter  Dose       Percentage 
Configuration (mr/hr-curie) (mr/hr-curie) Error 


Pi  0.0375  ±  0.0081      0.0264  ±  0.0008        42.0 

P2  0.0451  ±  0.0061      0.0291  ±  0.0009        55.0 

P3  0.0350  ±  0.0048      0.0240  ±  0.0007        45.8 


spectral  doses.   The  spectral  dose  is  equivalent  to  D  obtained  from  Eq.  (36) 
normalized  to  a  one  curie  source  strength.   The  per  cent  error  is  equal  to 
100%  (spectral  dose  -  dosimeter  dose) /dosimeter  dose. 

The  results  presented  in  Table  IV  indicate  an  average  percentage  differ- 
ence between  the  spectral  and  dosimeter  doses  of  47.6  per  cent.   The  percent- 
age error  on  the  individual  doses  ranged  from  13.5  to  21.6  per  cent  for  the 
spectral  doses  and  from  2.9  to  3.0  per  cent  for  the  dosimeter  doses.   These 
individual  errors  are  an  indication  of  the  precision  with  which  the  experi- 
mental measurements  were  made  and  the  error  introduced  in  reducing  the 
experimental  data. 

It  should  also  be  noted  that  the  difference  between  the  spectral  and 
dosimeter  doses  is  such  that  the  magnitudes  of  these  two  doses  are  not  within 
one  standard  deviation  of  each  other.   Therefore  it  is  apparent  that  there 
is  an  error  in  one  or  both  of  the  methods  used  to  determine  the  exposure 
doses  presented  in  Table  IV.   Since  the  range  of  the  percentage  errors  given 
in  Table  IV  is  not  large  it  is  not  unreasonable  to  suspect  that  the  error  is 
systematic. 
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At  this  point  it  must  be  recalled  that  one  of  the  primary  objectives  of 
this  work  was  to  determine  the  accuracy  of  the  unfolding  technique.   To 
achieve  this  goal  it  was  pointed  out  (see  Section  2.2)  that  first  the  accuracy 
of  the  dosimeter  results  must  be  established.   Thus  the  validity  of  the 
method  used  to  calculate  the  exposure  doses  from  the  dosimeter  data  must  be 
confirmed.   In  addition,  estimates  of  the  errors  associated  with  any  of  the 
approximations  used  in  determining  the  dosimeter  doses  should  be  obtained. 
If  it  can  be  sho^im  that  the  dosimeter  doses  presented  in  Table  IV  were 
obtained  by  a  valid  method  and  that  they  are  accurate,  within  some  maximum 
error  limit,  then  the  case  where  there  is  not  an  overlap  between  this  maximum 
error  limit  on  the  dosimeter  doses  and  the  error  limit  on  the  spectral  doses 
would  indicate  that  the  spectral  doses  were  in  error.   Therefore  an  estima- 
tion of  the  errors  associated  with  the  approximations  used  in  obtaining  the 
dosimeter  calibration  equations  as  well  as  a  detailed  discussion  concerning 
the  validity  of  these  calibration  equations  is  given  in  the  following  para- 
graphs . 

The  dosimeter  doses  presented  in  Table  IV  are  almost  entirely  dependent 
on  the  equations  for  the  calibration  curves  as  developed  in  Section  5.0. 
In  the  development  of  these  curves,  two  approximations  were  made:  (1)  that 
the  average  energy  of  the  gamma-ray  distribution  could  be  approximated  by 
the  single  scattered  energy  of  the  gamma  rays  scattered  from  the  concrete 
surface  and  (2)  that  the  curve  of  the  collimator  correction  factor  R  versus 
gamma-ray  energy  could  be  represented  by  a  straight  line.   Considering  the 
first  of  these  assumptions  it  was  previously  stated  that  the  single  scattered 
energy  would  be  higher  than  the  average  energy,  but  that  this  overestimation 
would  result  in  an  error  of  no  more  than  ten  per  cent.   The  second  assumption 
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could  result  in  an  error  of  approximately  five  per  cent  as  previously  stated 

in  Section  5.0.   Both  of  these  errors  together  could  then  possibly  increase 

the  confidence  limits  on  the  dosimeter  doses  by  as  much  as  15  per  cent. 

The  previous  discussion  would  tend  to  indicate  that  if  appropriate 

corrections  had  been  made  for  the  collimator  interactions,  as  pertains  to 

the  dosimeter  measurements,  that  the  dosimeter  and  spectral  doses  would  be 

in  closer  agreement.   However,  upon  closer  examination  of  the  equations  for 

the  calibration  curves  it  can  be  found  that  this  is  not  the  case.   The 

validity  of  this  statement  becomes  apparent  when  no  corrections  are  made  for 

the  collimator  interactions,  i.e.,  letting  (R  /R_  )  =  1.0  in  Eq.  (19).   If 

s  hi 

the  calibration  equations  of  Section  5.0  are  developed  under  this  condition 

the  average  percentage  difference  between  the  spectral  and  dosimeter  doses 

is  20  per  cent  and  the  spectral  doses  still  have  the  contribution  due  to 

collimator  interactions  removed.   Going  one  step  further,  if  one  calculates 

the  dosimeter  doses  presented  in  Table  IV  using  the  calibration  equations 

developed  in  Appendix  C  and   corrects   the  measured  dosimeter  responses  for 

the  angle  of  incidence  of  the  radiation  on  the  dosimeter,  the  doses  obtained 

are  within  four  per  cent  of  those  calculated  using  (R  /R^  )  =  1.0  and  this 

s   Ej 

difference  can  be  attributed  to  the  regression  analysis  used  in  developing 

these  two  calibration  equations. 

It  might  further  be  argued  that  the  fact  that  the  doses  obtained  using 

the  calibration  equations  developed  in  Section  5.0  converge  to  both  the  doses 

calculated  using  (R  /R   )  =1.0  and  using  the  calibration  equation  of  Appendix 

S   E  ]^ 

C  is  a  necessary  condition,  but  not  sufficient  to  substantiate  the  validity 
of  the  calibration  equations  of  Section  5.0.   Although  this  argument  is  valid 
in  that  the  doses  obtained  using  these  different  calibration  equations  must 
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all  converge  to  the  same  doses  (including  the  collimator  effects)  it  neglects 

one  Important  factor.   Since  both  calibration  curves  (that  of  Appendix  C  and 

those  of  Section  5.0)  converge  to  the  same  doses  vzhen  no  correction  is  made 

for  the  collimator  interactions  then  the  only  source  of  error  left  is  in  the 

ratio  R  /R^   or  in  the  source  strengths  of  the  cobalt-60  sources  used  in 
s   El  ^ 

obtaining  the  scattered  gamma-ray  distributions.   Now  the  ratio  R  /R   was 
determined  using  the  curve  of  the  collimator  correction  factor  R  versus  the 
gamma-ray  energy  described  in  Section  5.0  and  this  curve  is  identical,  within 
a  few  per  cent,  to  the  curve  used  in  making  corrections  for  the  collimator 
interactions  in  the  development  of  the  response  matrix.   Hence  a  change  in 
the  curve  of  the  collimator  correction  factor  versus  gamma-ray  energy  should 
result  in  an  approximately  equal  change  of  magnitude  in  both  the  spectral 
and  dosimeter  doses. 

With  reference  to  the  cobalt-60  source  used  in  the  dosimeter  measure- 
ments it  was  found  that  the  calibrated  value  of  82.6  +  1.8  curies  was  within 
statistical  agreement  of  the  value  quoted  by  the  vendor.   The  strength  of 
the  cobalt-60  source  used  in  the  spectral  measurements  (0.243  +  0.007  curies) 
was  approximately  twenty- two  per  cent  lower  than  the  vendor's  quoted  value. 
However  upon  calibrating  this  source  again  using  a  different  experimental 
configuration  and  detector,  a  value  for  the  source  strength  was  determined 
which  was  within  five  per  cent  of  the  value  obtained  in  the  previous  calibra- 
tion.  Thus  it  was  concluded  that  the  vendor's  value  for  this  source  was  in 
error  and  that  the  value  (0.243  +  0.007  curies)  used  in  this  report  was 
correct. 

There  is  one  further  point  which  must  be  considered  concerning  the 
applicability  of  the  calibration  curves  of  Section  5.0  to  the  dosimeter 
responses  obtained  from  the  scattered  gamma-ray  measurements.   In  particular, 
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consideration  must  be  given  to  the  difference  in  the  angular  distributions 
of  the  plane  source  used  in  developing  the  calibration  curves  and  the  simu- 
lated plane  sources  considered  in  the  scattered  gamma-ray  measurements.   Now 
the  response  of  the  dosimeter  is  proportional  to  the  magnitude,  as  well  as 
the  spatial  distribution,  of  the  gamma-ray  flux  throughout  the  volume  of  the 
dosimeter  and  these  quantities  are  in  turn  dependent  on  the  angular  distribu- 
tion of  the  plane  source.   The  magnitude  of  the  gamma-ray  flux  incident  on 
the  detector  is  taken  into  account  by  the  very  nature  of  the  calibration 
curves.   However  the  calibration  curves  do  not  apply  when  the  spatial  dis- 
tribution of  the  gamma-ray  flux  throughout  the  volume  of  the  dosimeter  is 
different  from  that  encountered  in  the  calibration  experiment  of  Section  5.0. 
Now  in  the  calibration  experiment  of  Section  5.0  a  plane  isotropic  source 
v/as  used  whereas  in  the  scattered  radiation  measurements  the  angular  distri- 
butions of  the  plane  sources  were  anisotropic  (17,18). 

Now  since  it  was  previously  stated  that  the  spatial  distribution  of  the 
gamma-ray  flux  was  dependent  on  the  angular  distribution  of  the  plane  source 
a  theoretical  investigation  of  this  problem  was  carried  out.   (See  Appendix 
E.)   This  investigation  was  designed  to  determine  if  the  differences  in  the 
angular  distributions  which  were  encountered  in  this  work  could  result  in  a 
significant  change  in  the  spatial  distribution  of  the  gamma-ray  flux  through- 
out the  volume  of  the  dosimeter.   Expressions  were  derived  for  the  spatially 
dependent  gamma-ray  flux  as  a  function  of  position  in  the  collimator  for 
various  plane  source  angular  distributions.   The  numerical  results  obtained 
from  this  investigation  indicated  that  for  the  angular  distribution  of  the 
scattered  radiation  sources  simulated  in  this  section  the  spatial  distribu- 
tion of  the  gamma-ray  flux  throughout  the  volume  of  the  dosimeter  was  nearly 
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the  same  as  that  resulting  from  the  isotropic  source  used  in  the  dosimeter 
calibration  experiment  of  Section  5.0.   Therefore  it  was  concluded  that  the 
difference  in  the  angular  distributions  of  the  plane  sources  used  in  the 
calibration  experiment  of  Section  5.0  and  those  encountered  in  this  section 
did  not  have  any  effect  on  the  applicability  of  the  calibration  curves  of 
Section  5.0  to  the  dosimeter  responses  obtained  in  the  scattered  gamma-ray 
measurements  of  this  section. 

From  the  preceding  discussion  it  appears  that  the  method  used  to  deter- 
mine the  dosimeter  doses  presented  in  Table  IV  was  valid.   In  addition,  it 
has  been  pointed  out  that  the  approximations  used  in  developing  the  equations 
for  the  calibration  curves  of  Section  5.0  could  result  in  a  total  error, 
the  magnitude  of  which  could  be  as  high  as  fifteen  per  cent.   Therefore  the 
maximum  error  associated  with  the  dosimeter  doses  presented  in  Table  IV  has 
a  range  of  from  17.9  to  18.0  per  cent.   Hence  the  previous  discussion  indi- 
cates that,  within  the  maximum  error  limits  presented  above,  the  dosimeter 
doses  tabulated  in  Table  IV  are  accurate  representations  of  the  exposure 
doses  resulting  from  the  experimental  configurations  considered  in  this  work. 

With  the  establishment  of  the  accuracy  of  the  dosimeter  doses  as  stated 
in  the  previous  paragraph  and  since  the  spectral  and  dosimeter  doses  are 
still  not  within  statistical  agreement  (their  error  limits  do  not  overlap), 
then  following  the  reasoning  developed  in  Section  2.2,  it  must  be  concluded 
that  there  is  an  error  in  the  spectral  doses  presented.   As  three  different 
experimental  configurations  were  considered  and  the  geometry  of  each  was 
reproduced  three  times  it  is  not  likely  that  such  an  error,  which  appears  to 
by  systematic,  is  due  to  the  experimental  spectral  measurements.   This 
indicates  that  the  error  is  in  the  unfolding  codes.   Previously  it  was 
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pointed  out  (see  page  42)  that  the  appearance  of  negative  values  in  the 
unfolded  pulse-height  distributions  as  illustrated  in  Figures  8,  9  and  10 
indicated  that  the  response  matrix  was  not  completely  valid.   Also  a  possible 
source  of  error  which  could  result  in  systematic  errors  such  as  evidenced  in 
this  work  is  the  incorrect  normalization  of  the  response  functions  used  in 
generating  the  response  matrix. 

Since  a  detailed  analysis  of  the  unfolding  codes  is  beyond  the  scope 
of  this  work  no  further  comments  can  be  made  concerning  the  source  of  errors 
in  the  unfolding  codes  on  the  basis  of  this  work.   However,  it  should  be 
noted  that  the  development  of  new  response  functions  and  a  new  response 
matrix  for  the  collimator  spectrometer  system  used  in  this  work  are  the 
subjects  of  two  research  projects  (19,20)  currently  underway  in  the  Depart- 
ment of  NucJ.ear  Engineering  at  Kansas  State  University. 
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7.0  CONCLUSIONS 

As  a  result  of  this  work  the  following  conclusions  can  be  made  concerning 
the  experimental  method  used  in  evaluating  the  unfolding  technique. 

1)  The  method  used  in  calibrating  the  dosimeter  for  polyenergetic 
radiation  using  a  calibrated  plane  source  of  cobalt-60  was  valid. 

2)  The  dosimeter  doses  obtained  using  the  calibration  referred  to 
in  1)  were  accurate  representations,  within  the  maximum  error 
limits,  of  the  true  exposure  doses  resulting  from  the  experimental 
configurations  considered. 

3)  The  experiment  designed  was  adequate  to  enable  a  determination  of 
the  accuracy  of  the  exposure  doses  computed  using  the  unfolding 
technique. 

Also  it  was  found  that  the  spectral  doses  obtained  using  the  unfolding 
codes  considered  in  this  work  were  larger  than  the  dosimeter  doses  and  were 
not  within  one  standard  deviation  of  the  dosimeter  doses.  Therefore,  since 
the  results  of  this  work  indicated  that  the  dosimeter  doses  were  accurate, 
within  their  maximum  error  limits,  it  was  concluded  that  the  spectral  doses 
were  in  error.  This  error  appeared  to  be  systematic.  Furthermore  it  was 
concluded  that  the  most  probable  source  of  this  error  was  in  the  normaliza- 
tion of  the  response  functions  and/or  the  generation  of  the  response  matrix 
itself. 

Although  a  detailed  analysis  of  the  unfolding  codes  was  beyond  the  scope 
of  this  work  the  results  of  this  work  do  indicate  that  such  an  analysis  of  the 
unfolding  codes  is  warranted.   At  the  present  tim.e  work  is  progressing  towards 
the  development  of  new  response  functions  and  a  new  response  matrix  (19,20) 
which  it  is  hoped  will  correct  for  the  shortcomings  of  the  response  matrix 
considered  in  this  work. 
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8.0   SUGGESTIONS  FOR  FURTHER  STUDY 

The  primary  purpose  of  this  work  was  to  develop  a  method  for  determin- 
ing the  accuracy  of  exposure  doses  calculated  using  an  unfolding  technique 
and  the  application  of  this  method  to  a  particular  unfolding  technique.   An 
obvious  extension  of  this  work  would  be  to  develop  an  experiment  to  determine 
the  accuracy  of  the  unfolded  pulse-height  distributions.   This  could  be 
accomplished  using  the  collimator-spectrometer  system  to  measure  scattered 
radiation  for  selected  gamma-ray  scattering  angles.   The  scattering  angles 
chosen  would  necessarily  be  those  for  which  Monte  Carlo  results  are  avail- 
able.  Then  the  unfolded  pulse-height  distributions  could  be  compared  against 
the  Monte  Carlo  results. 

The  results  of  this  work  clearly  indicated  that  the  unfolding  codes 
developed  for  the  collimator-spectrometer  system  used  in  this  work  should  be 
examined  in  detail.   Towards  this  end,  the  development  of  new  response  func- 
tions and  a  new  response  matrix  for  this  system  is  currently  in  progress  (19, 
20).   The  ability  of  the  unfolding  codes,  using  this  new  response  matrix,  to 
accurately  predict  exposure  doses  can  be  determined  using  the  experimental 
spectra  measured  for  this  work  and  comparing  the  resulting  spectral  doses  to 
the  dosimeter  doses  in  the  manner  outlined  in  this  report. 

The  dependence  of  the  collimator  correction  factor  on  the  gamma-ray 
energy  also  deserves  further  experimental  work.   As  was  previously  noted, 
Mather  has  investigated  this  problem  in  detail  on  a  theoretical  basis;  however, 
there  is  a  scarcity  of  experimental  results -in  the  literature  on  this  subject. 
The  relationship  between  the  collimator  correction  factor  and  the  gamma-ray 
energy  could  be  investigated  using  the  techniques  presented  in  this  work. 
This  would  involve  using  the  plane  simulator  with  a  number  of  monoenergetic 
sources  having  a  range  of  gamma-ray  energies. 
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Also  this  work  has  indicated  that  there  is  a  need  for  more  experimental 
work  on  the  angular-energy  dependence  of  the  Victoreen  ten  milliroentgen 
dosimeter.   An  experimental  study  of  this  nature  would  be  especially  valuable 
at  low  gamma-ray  energies.   However,  at  these  low  energies  one  of  the  great- 
est obstacles  is  the  acquisition  of  appropriate  monoenergetic  sources. 
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APPENDIX  A 
SPECTROMETER  QUALITY  CONTROL 

In  order  that  some  degree  of  confidence  could  be  placed  on  the  spectral 
data  obtained  using  the  gamma-ray  scintillation  collimator-spectrometer 
system,  a  quality  control  program  was  initiated  to  be  used  with  this  system. 
Although  quality  control  is  primarily  concerned  with  maintaining  the  quality 
of  a  system  at  a  certain  level  and  not  the  establishment  of  quality,  it  is 
apparent  that  through  the  correction  of  deficiencies  in  the  system  the 
quality  of  the  system  as  a  whole  may  be  improved. 

The  quality  control  of  scintillation  spectrometers  has  been  extensively 
investigated  (21,22,23,24)  and  the  results  have  indicated  that  drifts  in  the 
gain  and  baseline  may  completely  invalidate  the  calibration  of  the  system  and 
in  any  case  will  decrease  the  precision  of  any  measurements  which  might  be 
made.   These  drifts  result  in  fluctuations  in  the  energy  per  channel  and  the 
total  energy  range  being  examined  through  the  use  of  the  multichannel 
analyzer. 

It  is  of  prime  importance  to  maintain  a  constant  value  for  the  energy 
per  channel  and  the  total  energy  range  when  employing  a  scintillation  spec- 
trometer system.   However,  drift  in  the  baseline  (zero  energy  position) 
and/or  the  overall  gain  of  the  system  will  change  these  values.   Drift  due 
to  fluctuations  in  the  gain  is  primarily  attributed  to  small  changes  in 
(1)  the  conversion  and  optical  efficiency  of  the  crystal,  (2)  the  photomulti- 
plier  tube  gain,  and  (3)  the  linear  amplifier  gain.   Gain  drift  results  in  a 
change  of  the  energy  per  channel  of  the  system.   The  zero  energy  intercept 
or  baseline  is  determined  by  the  acceptance  threshold  of  the  m.ultichannel 
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analyzer.   Hence  a  baseline  drift  will  uniformly  shift  the  entire  spectrum 
as  well  as  relocate  the  position  of  the  zero  energy  channel. 

Fite  et  al.  (21)  have  developed  a  system  to  correct  for  the  drift  due 
to  gain  changes  and  changes  in  the  baseline  or  zero  point  of  the  analyzer 
system.   Their  method  utilizes  two  radioactive  sources  as  standards  and  the 
multichannel  analyzer  core  memory  is  then  monitored  to  determine  if  the 
peaks  are  in  the  proper  channels.   If  the  peaks  do  not  appear  in  the  proper 
channels  the  gain  and/or  baseline  are  altered  through  the  use  of  a  servo- 
system  to  bring  the  peaks  to  the  correct  channels. 

A  similar  system  has  been  developed  by  de  Waard  (22)  to  correct  for 
drifts  due  to  fluctuations  in  the  gain.   However,  he  considers  only  the 
drift  in  gain  and  indicates  that  perturbations  in  the  system  due  to  baseline 
drift  are  negligible  compared  to  those  caused  by  drift  in  the  gain. 
Corrections  for  the  drift  in  gain  are  made  through  the  use  of  a  feedback 
system  to  the  photomultiplier  tube  voltage. 

Covell  (23)  has  pointed  out  that  the  above  techniques  introduce 
extraneous  "background"  spectra  which  tend  to  complicate  the  interpretation 
of  the  measured  spectrum.   Therefore,  for  this  research,  a  control  procedure 
was  developed  which  utilized  a  daily  check  of  the  analyzer  system  using  a 
zinc-65  source  as  a  standard.   Corrections  were  then  made  manually  to  compen- 
sate for  any  changes  in  the  gain  and  baseline  demonstrated  by  the  spectrum 
for  the  standard  source.   A  similar  program  for  quality  control  has  been 
incorporated  by  Crouch  and  Heath  (24)  in  their  work.   Although  these  two 
methods  (23,24)  of  obtaining  a  control  over  the  analyzer  system  do  not  pro- 
vide a  continuous  monitor  of  the  system  they  do  eliminate  the  extraneous 
"background"  spectrum  noted  earlier. 
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In  order  for  the  quality  control  program  initiated  for  this  work  to  be 
effective  it  was  necessary  that  the  entire  procedure  (daily  sampling,  data 
evaluation,  logging  and  estimation  of  the  quality)  be  as  brief  as  possible 
and  still  maintain  a  high  level  of  efficiency. 

This  quality  control  program  is  in  agreement  with  Covell  (23)  and 
Crouch  and  Heath  (24)  insofar  as  drift  in  the  baseline  is  considered  to  be 
an  important  contributor  to  the  overall  drift  in  the  system  as  a  whole. 
Therefore  the  daily  check  on  the  quality  of  the  system  included  a  means  of 
correcting  the  drift  due  to  both  effects. 

The  standardization  of  the  spectrometer  system  was  repeated  as  often 
as  was  deemed  necessary  but  at  least  once  daily.   The  pulse  height  and  posi- 
tion of  each  photopeak  was  determined  from  the  spectrum  obtained.   The  posi- 
tions of  these  photopeaks  were  compared  with  the  preset  location  for  each 
peak.   If  drift  in  the  system  had  occurred,  the  gain  and/or  baseline  of  the 
multichannel  analyzer  were  altered  manually  to  return  the  peaks  to  the  proper 
channels. 

The  zinc-65  source  was  chosen  as  the  standard  since  it  emits  two 
readily  detected  and  well  separated  gamma  rays  (0.51  and  1.114  MeV)  in  the 
energy  region  of  interest  and  because  it  has  a  sufficiently  long  half- life 
(245  days).   Although  this  is  a  short  half-life  compared  to  standard  sources 
in  general,  it  was  long  enough  to  render  the  problem  of  source  decay  negli- 
gible over  the  period  of  its  use.   The  strength  of  the  source  and  its  posi- 
tion relative  to  the  detector  were  selected  such  that  there  would  be  no 
drift  induced  because  of  an  excessively  high  count  rate.   Care  was  taken  in 
positioning  the  source  in  order  that  a  constant  geometry  between  the  source 
and  the  scintillation  detector  might  be  maintained  from  day  to  day. 
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In  addition  to  the  daily  monitor  on  the  gain  and  baseline,  a  routine 
examination  was  made  daily  to  determine  if  the  multichannel  analyzer  were 
operating  correctly.   This  examination  consisted  of  several  tests  designed 
to  evaluate  whether  the  memory  location  and  computer  sections  of  the  analyzer 
were  performing  properly  and  was  carried  out  in  the  manner  outlined  in  the 
manufacturer's  operating  manual  (25). 

To  avoid  small  drifts  in  the  system  resulting  from  changes  in  the 
temperature  of  the  scintillation  detector  the  scintillation  crystal  and 
photomultiplier  tube  were  kept  at  approximately  106°F  by  employing  a  constant 
temperature  circulating  water  pump  to  force  heated  water  through  coils 
located  within  the  collimator  (Figure  2).   Because  of  the  large  mass  of  the 
collimator  (approximately  3500  pounds),  a  period  of  about  twenty-four  hours 
was  required  to  raise  the  temperature  of  the  collimator  to  the  desired  level. 
After  the  temperature  had  stabilized,  the  thermostatically-controlled  pump 
was  operated  continuously  to  maintain  the  temperatureof  the  system.   Although 
the  face  of  the  crystal  was  exposed  to  the  ambient  temperature  through  the 
two-inch  by  f ourteen-inch  collimator,  this  system  did  minimize  both  the  chance 
of  a  sudden  change  in  the  temperature  of  the  crystal-photomultiplier  tube 
assembly,  and  the  spectral  drift  which  would  accompany  such  a  temperature 
change. 

The  multichannel  analyzer  and  associated  instruments  were  kept  in  a 
closed,  air  conditioned  room.   This  was  done  to  decrease  the  possibility 
either  of  the  instruments  overheating  or  of  failure  of  the  readout  equip- 
ment because  of  dirt  in  the  mechanical  parts. 

The  quality  control  program  proved  to  be  of  great  aid  in  maintaining 
good  quality  control  for  the  gamma-ray  scintillation  spectrometer  system  as 
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compared  to  the  system  performance  previous  to  the  initiation  of  this  pro- 
gram.  Because  of  the  nature  of  the  experimental  measurements  made  with  this 
system  it  was  extremely  difficult  to  determine  if  the  system  was  operating 
properly  without  first  analyzing  the  experimental  data  obtained.   Therefore 
this  quality  control  program  v;as  particularly  valuable  in  that  it  provided  an 
early  indication  of  system  deficiencies  and  their  probable  location  in  the 
system,  thereby  saving  experimental  time  which  would  have  otherwise  resulted 
in  the  acquisition  of  highly  erratic  experimental  data. 
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APPENDIX  B 


CALIBRATION  OF  THE  COBALT-60 
SIMULATED  PLANE  ISOTROPIC  SOURCE 


The  plane  isotropic  source  of  cobalt-60  which  was  simulated  using  the 
line  source  and  the  plane  simulator  illustrated  in  Figures  3  and  4  was 
calibrated  in  the  configuration  of  Figure  11  using  the  Victoreen  ten- 
milliroentgen  full-scale  dosimeter  which  had  been  calibrated  in  the  manner 
described  in  Appendix  C.   Measurements  were  made  for  one  complete  simulation 
of  the  plane  source  and  two  source-to-detector  distances.   A  minimum  of 
three  measurements  were  made  for  each  of  the  two  distances.   From  an  inves- 
tigation of  the  background  for  a  time  equivalent  to  that  required  for  a 
simulation  of  the  source  plane,  it  was  concluded  that  the  dose  due  to  back- 
ground was  negligible.   The  dosimeter  was  read  using  the  Technical  Operations 
charger-reader  and  these  readings  in  microamperes  were  then  corrected  for 
temperature  and  pressure  using  Eq.  (C-1)  of  Appendix  C.   The  readings  were 
then  normalized  to  an  exposure  time  of  twenty-eight  minutes. 

The  experimental  dosimeter  responses  in  microamperes  were  averaged 
for  each  plane  source-to-detector  distance  and  converted  to  doses  in  milli- 
roentgens  using  Eq.  (C-19)  of  Appendix  C.   Experimental  dosimeter  responses 
in  microamperes,  corrected  for  temperature,  pressure  and  exposure  time,  and 
the  corresponding  doses  in  milliroentgens  for  each  of  the  two  cases  considered, 

are  presented  in  Table  B-1. 

* 

Using  the  doses  from  Table  B-1,  the  strength  of  the  plane  source  D   in 

milliroentgens  per  twenty-eight  minutes  was  determined  from 


^Although  Dp,  which  could  not  be  measured  experimentally,  has  the  units  of 
dose,  it  actually  represents  the  product  of  source  strength,  exposure  time 
and  the  flux  to  dose  conversion  factor.   This  unusual  terra  was  employed  for 
convenience  in  carrying  out  the  calculations. 
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D  = 
z 


^  d2A 


Area  of  plane  source 


(B-1) 


in  which  D  =  dose  in  mr  for  the  detector  at  a  distance  z  from  the  plane 
z  ^ 

source, 

d^A.  =  differential  area  on  the  plane  source,  and 

r  =  distance  f rom  d^A  to  the  detector  position  in  centimeters. 


Table  B-1.   Experimental  dosimeter  responses  and  their  corresponding 
average  doses  used  in  the  calibration  of  the  simulated 
plane  source. 


Plane  Source-to- 
Detector  Distance 
(feet) 


Dosimeter  Response 

Normalized  to  28  Minutes 

(microamperes) 


Average  Dose 
Normalized  to  28  Minutes 
(milliroentgens) 


12.0 


14.0 


25.50 
27.12 
25.70 
25.91 

19.35 
19.58 
19.57 


2.71  +  0.09 


2.02  +  0.07 


Assuming  that  the  plane  source  is  uniform  and  referring  to  Figure  11,  the 
preceding  equation  becomes 


D  = 


D 

Att 


rXi 


rX2 


dy 


dx 


(x^+y^+z^) 


+ 


dx 


(x^+y^+z^) 


(B-2) 


Integration  with  respect  to  x  yields 
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tan 
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+  tan 


-1 
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(y2+z2)l/2 


(B-3) 


and  solving  for  D  results  in 
P 


D   =  4itD 
P      z 


rY 


1      2 

dy 


2..2^1/2 


(y^+z^) 


tan 


-1 


Xi 


(y  2+^2)1/2 


+  tan 


-1 


X2 


,2^.2^1/2 


(y^+zO 


/  .  (B-4) 


The  integral  term  in  the  denominator  of  Eq.  (B-4)  was  evaluated  for  the 

two  values  of  plane  source-to-detector  distance  using  a  Gaussian  quadrature 

technique.   The  numerical  integration  was  programmed  for  the  IBM-1410  digital 

computer.   Substitution  of  the  appropriate  values  into  Eq.  (B-4)  for  each 

plane  source-to-detector  distance  resulted  in  a  value  of  D   for  each  of  the 

P 

two  cases . 

It  should  be  noted  that  Eq.  (B-4)  does  not  take  into  account  the  effects 
of  the  exponential  attenuation  in  air  or  the  buildup  factors  due  to  air  and 
ground  scattering.   Since  the  product  of  these  three  terms  could  be  con- 
sidered as  a  constant  value  over  the  area  of  the  plane  source,  an  average 
value  for  this  product  was  determined  (26,27)  for  each  of  the  plane  source- 
to-detector  distances.   The  values  used  were  1.053  and  1.039.   These  numbers 
correspond  to  source-to-detector  distances  of  fourteen  and  twelve  feet 
respectively.   Thus  the  values  of  D   obtained  from  Eq.  (B-4)  were  divided  by 
the  appropriate  value  (1.053  or  1.039)  to  correct  for  exponential  attenuation 
and  buildup  due  to  scattering.   These  two  values  were  then  averaged  to  yield 
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D  =  657  +  16  milliroentgens.  (B-5) 


This  value  of  D  is  proportional  to  the  source  strength  for  the  simulated 


plane  isotropic  source  of  cobalt-60. 
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APPENDIX  C 


CALIBRATION  OF  THE  DOSIMETER 
USING  A  POINT  SOURCE 


A  calibration  curve  was  developed  to  determine  the  relationship  between 
the  dosimeter  response,  i.e.,  the  microampere  reading  obtained  with  the 
Technical  Operations  charger-reader,  and  the  exposure  dose  in  milliroentgens 
when  the  dosimeter  was  exposed  to  a  point  source.   This  was  accomplished  by 
computing  a  regression  line  using  the  experimentally  measured  dosimeter 
responses  and  theoretically  calculated  exposure  dose. 

The  experimental  data  were  obtained  by  exposing  the  dosimeter  to  a 
calibrated  point  source  of  cobalt-60  for  various  time  increments.   A  calibra- 
tion board,  constructed  in  such  a  manner  as  to  minimize  gamma-ray  scattering, 
was  utilized  so  that  the  source-to-detector  distance  could  be  maintained  at 
a  constant  value.   This  calibration  board  was  suspended  above  ground  level  in 
order  to  decrease  ground  scattering.   This  experimental  configuration  is 
illustrated  in  Figure  12. 

The  data  were  corrected  to  a  temperature  and  pressure  of  22°C  and  760  mm 
of  mercury  by  multiplying  the  experimental  dosimeter  response  in  microamperes 
by  a  dimensionless  correction  factor  defined  by 


[273.0  +  0.555  (T-32.0)]  760.0  ,^^. 

"^  (295.0)  (25.40)  P  ^     ^ 


in  which  T  was  the  temperature  in  degrees  Fahrenheit  and  P  was  the  pressure 
in  inches  of  mercury  during  the  experimental  run. 

The  exposure  dose  rate  which  would  be  expected  from  the  experiment  can 
be  theoretically  computed  using 
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Figure  12.   Experimental  configuration  for  the  dosimeter 
calibration  using  a  point  source 
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in  which  P  Mt~)  =  theoretical  exposure  dose  rate  in  milliroentgens  per 

hour, 
K  =  conversion  factor  (mr/hr) / (mc/cm^) , 
S  =  strength  of  the  cobalt-60  source  in  millicuries, 
Bl    =   buildup  factor  for  air  scattering 
B2  =  buildup  factor  for  ground  scattering, 
y   =  total  gamma-ray  attenuation  coef  f  icient .  (cm"'^ )  for 

1.25  MeV  gamma  rays,  and 
X  =  source-to-detector  distance  in  centimeters. 

The  strength  of  the  source  on  the  date  of  the  experiment  is  calculated 
using 

S  =  S^  exp(-0.693  t/t^^^^  <^'^-3) 

in  which   S.  =  strength  of  the  cobalt-60  source  in  millicuries  on  the  date 
1 

of  its  calibration, 
t   =  time  in  years  between  the  date  corresponding  to  S,  and  that 
corresponding  to  S,  and 
t,,j   =   half-life  of  cobalt-60  in  years. 

Substitution  of  the  appropriate  values  into  Eq.  (C-3)  yields 

S  =  (6.67  mc)exp((-0. 693/5. 27)4.06)  =  3.91  mc  (C-4) 

where     S,  =  6.67  mc  on  October  31,  1961. 
1 

*The  units  for  this  conversion  factor  (mr/hr) / (7. 4x10  photons/cm^-sec) 
have  been  chosen  to  facilitate  evaluation  of  Eq.  (C-2)  with  S  expressed 
in  millicuries. 
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The  conversion  factor  K  was  based  on  the  definition  of  the  roentgen. 
One  roentgen  corresponds  to  the  dose  given  to  one  gram  of  air  when  83.8 
ergs  of  energy  are  absorbed  in  that  mass  of  air  (28).   Then  the  energy 
flux  per  roentgen  is  defined  by 

83.8  ergs 

(— )  .   r-cm'^ 
p  air 

in  which  (y  /p)  .   is  the  gamma-ray  energy  absorption  mass  attenuation 
coefficient.   Converting  from  ergs  to  MeV  yields 

5.23x10^  MeV 

( — )  .    r-cm'^ 
p   air 


Noxi7  let  (^1^  =   1.0  and  then  Eq .  (C-6)  can  be  written  as 
E 

5.23x10^  MeV 

1.0  r  =  -n- .  (C-7) 

( — )  .    cm^ 
p   air 


From  Grotenhuis  (29),  the  value  of  (y  /p)  .  for  1.25  MeV  gamma  rays  (this 
is  the  average  energy  of  the  two  gamma  rays  emitted  by  cobalt-60)  is  0.268 
cm  /gm.  Substituting  this  value  into  Eq.  (C-7)  and  converting  from  MeV  to 
curie-second  yields 

^  _      .5.23x10^  MeV,  ,    1  Y   i  rl  disintegration,  .1  curie   ,.  , 

^•°  ^  =  ^0.0268   cm^  ^  ^1.25  MeV^  ^ 2^ ^  ^JThOQ^^^ 

'  sec 

=  2.11  X  10-2  curie-sec  (^_3) 

cm'^ 

or  1  r  =  0.586  x  iq-S  curie-hr   ^  ^^_^^ 

cm^ 
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Thus  the  conversion  factor  K  is 

K  =  1.71  X  10^  (inr/hr)/(mc/cm2)  .  (C-10) 

The  source  and  detector  were  suspended  five  feet  two  inches  above 
ground  level  and  91.3  cm  apart  as  shown  in  Figure  12,   With  these  parameters 
and  using  the  experimental  curve  presented  by  Clark  and  Batter  (26) ,  the 
buildup  factor  for  ground  scattering,  B2,  was  found  to  be  1.00686. 

The  value  of  the  buildup  factor  due  to  air  scattering  was  obtained  from 
the  expression  (27) 

-o  1  n  _i_  n  m     0.0632  yx  .^  ^^. 

Bi  =  1.0  +  0.92  yx  e  .  (C-11) 

in  which  yx  is  the  source-to-detector  distance  in  mean  free  paths  of  air. 
Substitution  of  the  appropriate  value  of  yx  yields  B^  =  1.00552.   The  value 
for  the  exponential  attenuation  was 

g-yx  ^  Q^gg^  ^  (C-12) 

The  substitution  of  the  appropriate  values  into  Eq.  (C-2)  results  in  a 
theoretical  exposure  dose  rate  of 


V  El  f^\    -    (1»  71x10  ^)  (3.91)  (1.00552)  (1.00686)  (0.994)  _  .  , ,  mr    ,   „. 
Th'^hr^  4tt(91.3)^  ^'        hr  *  ^        "^^ 


The  experimental  dosimeter  responses  and  the  theoretical  doses  are  presented 
in  Table  C-1.  The  experimental  dosimeter  responses  listed  in  this  table  have 
been  corrected  for  temperature  and  pressure  according  to  Eq.  (C-1). 

It  is  imperative  at  this  point  to  have  a  basic  understanding  of  the 
dosimeter  and  the  charger-reader  if  a  curve  is  to  be  fitted  to  the  data 
presented  in  Table  C-1.  In  the  use  of  this  detection  system  the  dosimeter  is 
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Table  C-1.   Data  used  in  the  calibration  of  the 
dosimeter  using  a  point  source 


Exposure  Time         Dosimeter  Response         Theoretical  Dose 
(minutes) (microam.peres  ) (milliroentgens) 


15.0  15.53  1.60 

15.53 
16.05 

25.0  25.84  2.67 

25.27 
26.31 
25.81 

35.0  36.16  3.74 

35.65 
35.62 


initially  charged  to  a  given  voltage  corresponding  to  a  zero  reading  on  the 
charger-reader  meter  scale.   When  the  dosimeter  is  exposed  to  ionizing 
radiation  electrons  are  ejected  from  the  dosimeter  walls  which  in  turn 
result  in  the  production  of  ion  pairs  in  the  dosimeter  chamber  gas.   These 
ion  pairs  are  then  sv;ept  away  to  the  detector  electrodes  due  to  the  poten- 
tial gradient  across  the  dosimeter.   Upon  reaching  the  electrodes,  the  ion 
pairs  have  the  effect  of  neutralizing  a  portion  of  the  original  dosimeter 
charge.   When  the  dosimeter  is  read  using  the  charger-reader  the  output 
voltage  is  the  difference  between  the  dosimeter  voltage  and  the  original 
charging  voltage.   This  output  voltage  is  thus  a  linear  function  of  the 
charge  collected,  i.e.,  the  integrated  dose  received  by  the  dosimeter. 

In  the  interest  of  clarity  and  with  the  preceding  discussion  in  mind, 
a  word  of  explanation  is  in  order  concerning  the  units  on  the  charger-reader 
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meter  scale.   The  units  of  this  scale  and  thus  the  units  associated  with  the 
dosimeter  responses  presented  in  this  report  are  microamperes.   However,  in 
the  preceding  paragraph  it  was  stated  that  the  quantity  measured  was  a 
voltage.   Thus  there  appears  to  be  a  contradiction  between  the  quantity- 
measured  and  the  quantity  which  is  obtained  from  the  charger-reader  meter 
scale.   This  apparent  contradiction  results  from  the  fact  that  when  the  manu- 
facturer assembled  the  charger-reader  an  ammeter  was  modified  in  order  to 
make  it  voltage  rather  than  charge  sensitive  but  the  appropriate  change  was 
not  made  in  the  units  on  the  scale  of  this  meter.   However,  because  the 
dosimeter  and  charger-reader  must  be  calibrated  before  being  used  and  because 
the  needle  deflection  on  the  charger-reader  meter  scale  is  a  linear  function 
of  the  dose  received  by  the  dosimeter  the  units  given  on  this  scale  are  not 
really  significant.   Thus  in  order  to  avoid  confusing  the  reader  who  is 
familiar  with  this  particular  charger-reader  the  author  has  elected  to  report 
the  dosimeter  response  in  terms  of  the  units  given  on  the  meter  scale,  i.e., 
microamperes. 

From  the  preceding  discussion  concerning  the  linearity  between  the 
deflection  of  the  needle  on  the  charger-reader  meter  scale  and  the  integrated 
exposure  dose  received  by  the  dosimeter,  then  if  the  charger-reader  meter  is 
zeroed  preceding  the  charging  of  the  dosimeter,  a  zero  exposure  of  the  dosim- 
eter to  ionizing  radiation  should  result  in  a  zero  microampere  reading  on 
the  charger-reader.   Following  a  check  of  the  charger-reader  and  dosimeter 
it  was  observed  that  for  no  exposure  to  radiation  the  dosimeter  response  was 
indeed  zero.   Also  upon  investigation,  it  was  learned  that  the  drift,  i.e., 
loss  of  initial  charge  due  to  leakage  and  other  effects,  in  the  dosimeter 
charge  was  less  than  0.5  microamperes  per  hour.   Thus  it  was  concluded  that 
loss  of  charge  because  of  drift  was  negligible. 
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Therefore  with  the  knowledge  that  the  response  of  the  dosimeter  to 
ionizing  radiation  is  a  linear  function  of  the  integrated  dose  to  which  the 
dosimeter  was  exposed  and  that  for  no  exposure  to  radiation  the  dosimeter 
response  is  zero,  the  data  of  Table  C-1  were  fitted  to  a  straight  line 
passing  through  the  origin  using  the  regression  analysis  outlined  in  Appendix 
D.   The  experimental  dosimeter  responses  were  taken  as  the  y  values  and  the 
theoretical  doses  as  the  x  values  in  this  analysis.   The  values  of  x  had  a 
constant  percentage  error  equivalent  to  that  associated  with  the  calibration 
source  strength.   Thus  the  Berkson  model  (30)  for  a  controlled  experiment  was 
utilized.   (See  Appendix  D.)   From  Eq.  (D-6)  of  Appendix  D  the  slope  of  the 
line  is 


(C-14) 


10 

I 

i=l 

x.y. 

10 

I 

i=l 

(x.)^ 

The  estimate  of  the  variance,  a  ,  from  Eq.  (D-20)  of  Appendix  D  is 

(s')2  =  s2  +  T^Z  v[x]  (C-15) 


in  which 


10 
(  I     x.y.)2 

1=1 


s'z  ^  ±:i± =  0.148  (C-16) 

£  -  1 


,2  -  i=l 


and 


b2v[x]  =  b2[0.05(^)]''  =  (0.05  y')2  (C-17) 
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where  y'  is  the  mean  of  m   new  observations  and  the  0.05  represents  the 
standard  deviation  in  the  calibration  source  strength. 

Since  the  calibration  curve  was  to  be  used  to  convert  from  microamperes 
to  milliroentgens ,  the  least  squares  line  in  reverse  was  determined.   To 
calculate  the  value  of  x  corresponding  to  some  y',  the  equation 

v'    ts'   1      (v'")2     1/2 

i=l 

can  be  used  if  the  condition  on  g' ,  as  defined  by  Eq.  (D-22)  of  Appendix  D, 
is  satisfied.   The  value  of  t  is  1.05  for  68%  confidence  (one  standard  devia- 
tion) and  nine  degrees  of  freedom.   The  value  of  g'  in  the  range  of  interest 

is  much  less  than  0.1  and  thus  the  least  squares  line  in  reverse,  as  given 

10 

above,  is  valid.   Substitution  of  the  values  of  b,  s',  t  and  I      (x.)^ 

i=l   ^ 

as  calculated  from  the  data  given  in  Table  C-1  yields 

x  =  0.104  y'  +  0.109[0.148  +  (0.05  y')2]-^^^[-  +  0.000138(y '  )2]  ^/^.  (C-19) 
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APPENDIX  D 
REGRESSION  ANALYSIS 

The  method  of  least  squares  is  a  mathematical  technique  for  fitting  an 
"optimum"  curve  through  a  given  number  of  points.   This  optimum  curve  is 
obtained  by  determining  the  parameters  of  the  curve  in  such  a  way  as  to 
minimize  the  sum  of  the  squares  of  the  deviations  between  the  values  of  the 
observed  and  predicted  points.   In  fitting  a  straight  line  which  passes 
through  the  origin  to  a  set  of  data  points,  the  following  two  cases  must 
be  considered;  1)  one  variable  subject  to  error  and  the  other  with  negli- 
gible error  and  2)  both  variables  subject  to  error. 

Following  the  method  outlined  by  Brownlee  (31) ,  the  case  where  y  is 
observed  with  error  and  x  is  assumed  to  be  free  of  error  is  considered  first. 
In  addition,  it  is  assumed  that  the  observed  values  of  y  have  a  normal  dis- 
tribution about  the  expected  value  y  with  a  variance  a   and  that  all  obser- 
vations are  independent. 

If  the  true  equation  of  the  line  is  represented  by 

Y  =  6x  (D-1) 

then  the  least  squares  analysis  of  the  data  will  yield  estimates  b  and  s^ 
of  S  and  o^.      Also  this  analysis  will  provide  a  method  of  determining  the 
distribution  of  these  estimates.   Assume  that,  for  each  value  of  x.,  i=l, 
2,  ...,  k, there  are  n.  values  of  y..   The  data  points  can  be  considered  by 
pairing  each  value  of  y.  with  its  corresponding  value  of  x..   For  example, 
if  k=2,  ni=2,  and  n2=3,  then  the  points 

(yil.xi)(yi2,xi)(y2i,X2)(y22.X2)(y2  3>X2)  (D-2) 
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can  be  represented  as 

(yi.xi)(y2,X2)(y3,X3) (y4,X4)(y5,X5)  ,  (D-3) 

where     x^  =  X2  and  X3  =  x^  =  X5. 

Now  if  the  estimated  least  squares  line  is  represented  by 

Y  =  bx  , 

then  the  sum  of  the  squares  of  the  deviations  between  the  observed  values 
y.  and  the  estimated  values  Y.  is 

£  a 

Q  =   I   (y,  -  Y  )2  =  I       (y   -  bx  )2  (D-4) 

1=1  1=1 

in  which  £  is  the  total  number  of  observations.   The  minimum  of  Q  is  obtained 
by  taking  the  first  derivative  of  Q  with  respect  to  b  and  equating  it  to  zero, 


i.e.,  ^=  -^     I      ^^i  -  ^^i>^i  =  0'  (°-^ 

i=l 


y   x.y, 
from  which  b  =  — .  (D-6) 


i=l 


The  estimate  of  a^  is  defined  as 


,2 


s2  =  ^ ,  (D-7) 


a  -  1) 

where  £  -  1  represents  the  degrees  of  freedom.   Now 

y.  =  (y.  -  Y.)  +  Y.  =  (y.  -  Y.)  +  bx.  .  (D-8) 

1     1    1     1     1    1      1 
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Squaring  this  equation  and  summing  over  i  yields 


Thus 


I      (7^)2  =1   (Yi  -  Y.)2  +  -f^ .  (D-9) 

i=l 


I      (y,)'  -[(  I   X  y  )2/  I      (x  )2] 

32  =  izl izl izl .  (D_10) 

£  -  1 


The  least  squares  line  can  be  used  in  reverse  to  calculate  a  value  of  x 
corresponding  to  some  y',  the  mean  of  m  new  observations.   The  point  estimate 
of  X  is  obtained  from 

X  =  ^'  .  (D-11) 

The  value  of  y'  which  would  correspond  to  this  value  of  x,  using  the  true 
least  squares  line,  is  n. 

From  Eq.  (D-1) ,  the  value  of  x  corresponding  to  this  value  of  n  is 


n 

X  =  — 


(D-12) 


and  if  this  true  value  of  x  is  designated  by  a,  then 

n  -  a6  =  0  .  (D-13) 

Define  the  variable  w  as 

w  E  y'  -  ba  .  (D-14) 

The  expected  value  of  w  can  be  represented  by 

E[w]  =  n  -  aE[b]  =  n  -  a6  =  0  (D-15) 
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and  its  variance  as 


V[w]  =  V[y']  +  a'^V[b]  =  -  + 


2        2  2 


m     £ 
i=l 


(D-16) 


Using  the  estimate  of  o^      given  by  Eq.  (D-10) ,  the  value  of  x  corre- 
sponding to  y'  is 


by' 


ts 


X  = 


b^  - 


T^^ 


T^^ 


I  (\y 

i=l 


i=l 


b^  - 


t2s2 


i=l 


_     1^   /y') 

m    £ 
i=l 


1/2 


(D-17) 


where  t  is  the  fractional  point  of  the  t  distribution  based  on  (£-1)  degrees 
of  freedom.   The  preceding  equation  for  x  is  governed  by  a  term  g  which  is 
defined  as 


t2s2 


i=l 


(D-18) 


If  the  value  of  g  is  less  than  0.1,  then  the  approximation 


b2  - 


t2s2 


I      (Xi)2 
i=l 


1  - 


t2s2 


b2  I       (X.) 

i=l 


2 


b2(l-g)  -  b' 


can  be  used  in  Eq.  (D-17).   Thus,  for  the  case  v;here  g  is  less  than  0.1,  the 
value  of  X  corresponding  to  y'  can  be  determined  from 


y   ,  ts 
^  =  b  ^"b 


1+     , 

m       £ 


iyJL) 


m2 


1/2 


(D-19) 


b2  I       (x)2j 
i=l    ^ 
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The  case  where  both  variables  are  subject  to  error  can  be  treated  by 
using  the  Berkson  model  for  controlled  experiments  (30,31).   Berkson  noted 
that  when  an  experiment  was  carried  out  in  such  a  manner  that  the  x  variable 
was  a  "controlled"  quantity  (i.e.,  for  each  observation  of  y.,  the  corre- 
sponding value  of  X  was  assigned  a  value  x.,  which  was  as  close  to  the  true 
value  of  X,  as  was  experimentally  possible),  then  the  least  squares  line 
could  be  determined  in  the  same  manner  as  was  previously  described  (one 
variable  subject  to  error).   However,   for  this  case,  the  variance  associated 
with  the  least  squares  line  must  contain  a  term  which  takes  into  account  the 
error  associated  with  the  x  variable. 

For  a  least  squares  line  through  the  origin,  the  slope  of  the  line  can 
be  found  using  Eq.  (D-6) .   The  estimate  of  a^  is  represented  by  Mandel  as 

(s')2  =  s2  +  b2v[x]  (D-20) 

where  s   is  defined  by  Eq.  (D-ID)  and  V[x]  is  the  variance  associated  with 
the  "controlled"  variable  x. 

In  order  to  use  the  least  squares  line  in  reverse  to  find  a  value  of  x 
corresponding  to  the  mean  of  m  new  observations,  y',  Eq.  (D-17)  should  be 
used.   However,  the  value  of  s  in  Eq.  (D-17)  must  be  replaced  by  the  s'  in 
Eq.  (D-20).   The  expression  for  x  is  then 


v'  ^  ts' 
X  =  -^  + 


b  -  b 


n 


1  I   (y')' 

m       ^ 

h'      I      (x)2J 
i=l 


1/2 

(D-21) 


under  the  condition  that  g'  is  less  that  0.1.   The  value  of  g*  is  determined 
from 

g.  =  __t_2(^:Jl_  .  (D-22) 


i=i    ^ 
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APPENDIX  E 


UNCOLLIDED  GAM>IA-RAY  FLUX  IN  THE  COLLIMATOR 
RESULTING  FROM  VARIOUS  PLANE  SOURCE  ANGULAR  DISTRIBUTIONS 


This  appendix  is  devoted  to  a  theoretical  investigation  of  the  varia- 
tion in  the  uncollided  gamma-ray  flux  within  a  collimator  resulting  from 
changes  in  the  angular  distribution  of  a  plane  source  located  at  one  end  of 
the  collimator.   In  particular  this  study  was  intended  to  determine  if  the 
spatial  distribution  of  the  flux  throughout  the  collimator  varied  as  the 
angular  distribution  of  the  plane  source  was  changed,  and  if  so,  the  extent 
of  this  variation.   The  derivation  of  formulae  for  the  uncollided  flux  as  a 
function  of  position  within  the  collimator  was  carried  out  using  a  general 
collimator-plane  source  geometry.   This  collimator-plane  source  geometry  is 
illustrated  in  Figure  13  and  consists  of  a  right  circular  cylindrical  colli- 
mator with  a  plane  source  positioned  so  that  it  is  perpendicular  to  the 
collimator  axis  and  adjacent  to  the  base  of  the  collimator. 

The  uncollided  gamma- ray  flux,  as  a  function  of  position  in  the  colli- 
mator, can  be  represented  by  an  integral  of  the  source  angular  distribution 
over  the  solid  angle  defined  by  the  collimator,  i.e., 


(R,z,(j))  = 


S(rj)  sec  e  d^^  (E-1) 


where  S(^)    is  the  source  angular  distribution,  d^.Q  is  the  differential  solid 
angle  and  the  angle  6  is  defined  in  Figure  13.   Since  only  azimuthally 
independent  source  functions  were  considered  then  3(9.)    =   S(9)  and  hence  the 
flux  was  also  azimuthally  independent  so  <|)(R,z,(j))  =  <I>(R,z).   It  should  be 
noted  that  Eq.  (E-1)  neglects  both  exponential  attenuation  and  buildup  due 
to  air  scattering. 
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Source  Plane 


Collimator  Wall 


Figure  13.   Collimator  source  geometry  used  in  the  derivation  of  the 
uncollided  flux  as  a  function  of  position  P  within  the 
collimator 
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The  differential  solid  angle  can  be  written  as 


d^^   =  — 7  cos 
r^ 


p  dp  dd) 

T— ^  COS 


(E-2) 


where  d^A,  r,  p,  <})  and  6  are  defined  in  Figure  13.   Also  from  Figure  13 
the  following  relation  can  be  obtained 


r^  =  z^  +  r2  +  p2  -  2Rp  cos  (^    . 


(E-3) 


Substitution   of   Eq.     (E-2)    into   Eq.     (E-1)    and   using   Eq.     (E-3)   yields 


<^(R,z)    = 


pdp 


27T 


S(9)    d(}) 

[z^+R^+p^   -   2Rp    cos   4)] 


(E-4) 


The  three  source  functions  considered  were 


s^e) 

1 

47T 

s"(e) 

3 

47T 

2 
COS 

S"^6) 

5 

4jr 

COS^ 

(E-5a) 
(E-5b) 
(E-5c) 


where  the  angle  6  is  defined  in  Figure  13.  These  source  functions  are 
normalized  such  that  integration  of  these  functions  over  47t  steradians 
yields  unity.  Eq.  (E-5a)  represents  an  isotropic  distribution  whereas 
Eqs.  (E-5b)  and  (E-5c)  represent  distributions  which  are  peaked  in  the  forward 
direction, i. e. ,  the  direction  perpendicular  to  the  source  plane.  From 
Figure  13  the  following  relation  can  be  obtained, 


cos  0  = 


[z2+R2+p2  _  2Rp  cos  ())] 


1/2 


(E-6) 
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The  fluxes  resulting  from  the  three  source  distributions  given  by  Eqs. 
(E-5a),  (E-5b)  and  (E-5c)  are  denoted  by  ^^(R,z),    $''"■'"  (R,z)  and  ^^^^iR,z) 
respectively.   Substitution  of  these  source  distributions  into  Eq.  (E-4) 
and  using  the  relationship  given  by  Eq.  (E-6)  yields 


fR 


^    (R,z) 


1_ 

Ait 


pdp 


2^ 


d(^ 


[z^+R^+p-^    -    2Rp    cos    (J)] 


(E-7a) 


4TT 


R 


pdp 


2-n 


d(j) 


[z^+R^+p^   -    2Rp    cos    <t>] 


(E-7b) 


R 


pdp 


2tt 


d(^ 


[z^+R^+p^   -    2Rp    cos    (J)] 


(E-7c) 


Eqs.    (E-7a) ,    (E-7b)    and   E-7c)   were  made  dimensionless   through   the   transforma- 
tion of  variables   given  by 


R  _        z  p 

^   =  R      '      ^   =  R      '      '^   =  R 
o  o  o 


(E-8) 


Application  of   this    transformation  of  variables    to   the   flux  equations  yields 


*^v,o     -i^ 


ydy 


2lT 


d()) 


[C   +v'^+lJ^   -    2viJ    cos    cf)] 


(E-9a) 


ydy 


27T 


d(() 


[^'^H-v^+y^  -   2vy   cos    tf)]' 


(E-9b) 


4)         (v,ri) 


5i' 
Att 


r2TT 


ydy 


d(|) 

[5^+v^+y^   -   2vy   cos   (fi]^ 


(E-9c) 
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Carrying  out  the  integrations  in  the  expressions  for  the  fluxes  given  above 
yields 

*  (v,ii)    =  ^  In  t 2F2 1  (E-lOa) 


^^  [(^2+1)2  +  v^  +  2v2(e2-l)]l/2 

(^2  +  ^Zp  I  (E-lOc) 


From  these  three  expressions  the  uncollided  gamma-ray  flux  can  be  determined 
for  any  point  in  the  collimator  in  terms  of  the  dimensionless  parameters  v 
and  E,   v;hich  are  directly  proportional  to  the  radial  and  axial  position  of  the 
point  in  the  collimator. 

Now  consider  the  collimator-plane  source  geometry  used  in  the  dosimeter 
calibration  of  Section  5.0  and  that  used  for  the  scattered  gamma-ray  measure- 
ments of  Section  6.0.   In  these  two  cases  the  source  plane  was  not  located 
adjacent  to  the  base  of  the  collimator,  however,  as  Eq.  (E-1)  was  based  on 
an  integration  over  the  solid  angle  the  distance  from  the  source  to  the 
point  of  interest  is  not  a  factor  in  determining  the  flux  as  long  as  the 
solid  angle  remains  constant. 

In  the  dosimeter  calibration  of  Section  5.0  a  plane  isotropic  source  was 
used.   On  the  other  hand,  theoretical  studies  (17,18)  indicate  that  the 
angular  distribution  for  radiation  scattered  from  a  concrete  slab  is  not 
isotropic.   In  fact  for  the  experimental  configurations  considered  in  this 
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work  it  appears  that  the  angular  distribution  of  the  scattered  radiation  is 
proportional  to  sec  6  where  6  is  the  angle  that  the  scattered  radiation  makes 
with  the  perpendicular  to  the  scattering  plane.   Now  recall  that  for  the 
scattered  radiation  measurements  the  scattering  surface  was  not  perpendicu- 
lar to  the  collimator  axis  but  rather  it  was  at  an  angle  of  72°  to  the 
collimator  axis  in  each  of  the  experimental  configurations  considered. 
Then  based  on  this  fact  and  the  previously  mentioned  theoretical  studies 
found  in  the  literature,  the  effective  source  seen  by  a  detector  located  in 
the  collimator  should  have  an  angular  distribution  which  is  peaked  in  the 
direction  of  the  detector.   Although  the  analytic  form  of  this  angular 
distribution  is  not  known,  the  theoretical  studies  (17,18)  indicate  that  it 
is  not  as  highly  peaked  in  the  direction  of  the  detector  as  the  cos  9 
distribution  considered  in  the  derivation  of  Eq.  (E-lOc).   Hence  it  appears 
that  Eq.  (E-lOa)  is  applicable  to  the  calibration  experiment  and  that 
Eq.  (10-c)  should  provide  an  approximate  upper  limit  with  respect  to  the 
peaking  of  the  angular  distributions  encountered  in  the  scattered  gamma-ray 
measurements. 

Now  since  the  purpose  of  this  appendix  was  to  determine  if  the  spatial 
distribution  of  the  uncollided  gamma-ray  flux  varied  as  the  angular  distri- 
bution of  the  plane  source  was  changed,  Eqs .  (E-lOa) ,  (E-lOb)  and  (E-lOc) 
were  evaluated  for  different  positions  in  the  collimator.   In  order  that  the 
results  of  these  computations  could  be  used  in  support  of  the  experimental 
techniques  used  in  this  work,  the  flux  was  determined  for  various  points  in 
a  particular  right  circular  cylindrical  volume  of  the  collimator.   This 
volume  corresponds  to  the  volume  occupied  by  the  dosimeter  during  the  experi- 
mental measurements  of  Sections  5.0  and  6.0.   (See  Figure  6.)   Calculations 


93 


were  made  for  points  along  the  centerline  (v=0) ,  edge  (v=l)  and  base  of  this 

cylindrical  volume.   So  that  the  results  from  the  different  plane  source 

angular  distributions  could  be  readily  compared,  the  fluxes  at  the  various 

points  considered  were  normalized  to  the  flux  at  the  center  of  the  volume 

$   for  each  of  the  source  distributions.   Table  E-I  is  a  tabulation  of  these 
c 

normalized  fluxes  for  the  three  plane  source  angular  distributions  considered. 


Table  E-I.   Normalized  fluxes  as  a  function  of  position  in  the 
collimator  and  plane  source  angular  distribution. 

$  =  =  0.86356  X  10~   gammas/cm^-sec 


II 


III 


=  2.5862  X  10  ^  gammas /cm^-sec 


$    =  4.3029   X  10  ^  gamraas/cm^-sec 


■^  (v=0,C) 


c, 

c 

15.625 

1.1834 

16.175 

1.1044 

16.725 

1.0331 

17.275 

.96847 

17.825 

.90972 

18.375 

.85615 

C 

$-^(v=l,0 
c 

15.625 

1.1786 

16.175 

1.1002 

16.725 

1.0294 

17.275 

.96525 

17.825 

.90688 

18.375 

.85364 

c 


c 


1.1830 
1.1042 
1.0330 
.96852 
.90986 
.85636 


^^\v-l,g) 
c 


1.1734 
1.0959 
1.0257 
.96211 
.90419 
.85134 


1.1826 
1.1040 
1.0330 
.96857 
.91001 
.85658 


III 


(v=l,C) 


$ 


III 


1.1684 
1.0916 
1.0221 
.95898 
.90153 
.84907 
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Table  E-1  (continued) 


0.0 
0.2 
0.4 
0.6 
0.8 
1.0 


I  II  ITT 

(v,C=15.625)        J)   (v,C=15.625)        'I'   (v,C=15.625) 


c 

c 

c 

1.1834 

1.1830 

1.1826 

1.1832 

1.1826 

1.1820 

1.1826 

1.1815 

1.1803 

1.1816 

1.1795 

1.1775 

1.1803 

1.1769 

1.1735 

1.1786 

1.1734 

1.1684 

From  a  comparison  of  the  values  of  <J>  ,  $   ,  and  $    given  in  Table  E-1, 
it  can  be  seen  that  the  flux  increases  as  the  angular  distribution  of  the  plane 

source  becomes  more  highly  peaked  in  the  direction  perpendicular  to  the 

I       I 
source  plane.   However,  the  values  of  the  normalized  fluxes  $  (v,C)/$  , 

*   ('v^>?.)/$   and  $    (v,^)/^    for  any  point  in  that  portion  of  the  collima- 
tor volume  under  consideration  were  in  close  agreement  to  one  another.   In 
fact  a  comparison  of  these  three  ratios  for  any  of  the  points  (v,0  considered 
in  Table  E-1  alxv^ays  resulted  in  an  agreement  between  them  which  was  less  than 
one  per  cent.   Therefore  it  can  be  concluded  that  for  the  collimator-plane 
source  systems  considered  in  this  appendix  a  change  in  the  angular  distribu- 
tion of  the  plane  source  results  in  a  negligible  change  in  the  spatial  distri- 
bution of  the  flux  within  the  collimator  volume  considered. 

Now  it  has  been  pointed  out  that  the  plane  isotropic  source  case  was 
applicable  to  the  dosimeter  calibration  of  Section  5.0  and  that  the  distri- 
bution should  provide  an  upper  limit  with  reference  to  the  peaked  angular 
distributions  of  the  effective  plane  sources  in  the  scattered  radiation 
measurements  of  Section  6.0.   Also  the  numerical  results  presented  in 
Table  E-1  indicated  that  the  spatial  distribution  of  the  flux  throughout  the 
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volume  which  the  dosimeter  would  occupy  is  the  same  for  the  three  angular 
distributions  considered.   Thus  the  conclusions  resulting  from  this  study 
can  be  extended  to  the  work  done  in  the  remainder  of  this  report  by  stating 
that  in  both  the  dosimeter  calibration  of  Section  5.0  and  the  scattered 
radiation  measurements  of  Section  6.0  the  spatial  distribution  of  the  un- 
collided  gamma-ray  flux  throughout  the  volume  of  the  dosimeter  was  essentially 
the  same.   Hence  it  can  be  concluded  that  the  application  of  the  calibration 
curves  developed  in  Section  5.0  and  used  to  analyze  the  dosimeter  responses 
obtained  from  the  scattered  radiation  measurements  of  Section  6.0  was  valid 
even  though  the  angular  distributions  of  the  plane  sources  considered  in 
Sections  5.0  and  6.0  were  not  the  same. 
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ABSTRACT 

A  study  was  made  to  determine  experimentally  the  validity  of  using  an 
unfolding  technique,  developed  for  the  analysis  of  data  obtained  using  a 
collimator-spectrometer  system,  in  the  computation  of  exposure  doses.   This 
determination  was  based  on  a  comparison  between  exposure  doses  computed  from 
experimental  complex  gamma-ray  energy  spectra  using  the  unfolding  technique 
and  exposure  doses  resulting  from  measurements  with  an  air  equivalent  ioniza- 
tion chamber  (dosimeter). 

The  experimental  phase  of  this  work  consisted  of  measuring  exposure 
doses  for  various  gamma-ray  energy  distributions  using  both  an  Nal(Tl)  scin- 
tillation detector  and  a  dosimeter.   An  energy  distribution  of  gamma  rays  was 
obtained  by  scattering  cobalt-60  radiation  from  a  concrete  slab.   The  two 
detection  instruments  were  placed  (alternately)  in  a  lead  collimator  to 
enable  the  measurement  of  gamma-ray  energy  distributions  corresponding  to 
three  different  experimental  configurations. 

An  expression  was  derived  relating  the  response  of  the  dosimeter  to  a 
plane  monoenergetic  source  to  its  response  to  a  plane  polyenergetic  source 
of  equivalent  strength  and  known  energy  distribution.   Using  this  expression 
calibration  equations  relating  dosimeter  response  to  exposure  dose  were 
obtained  for  each  of  the  energy  distributions  to  be  considered  by  exposing 
the  dosimeter  to  a  simulated  plane  isotropic  source  of  cobalt-60.   The 
validity  of  this  calibration  technique  for  the  dosimeter  was  examined  in 
detail. 

As  a  result  of  the  comparison  between  the  exposure  doses  obtained  using 
the  dosimeter  and  those  obtained  by  applying  the  unfolding  technique  to  the 


complex  spectra  measured  it  v/as  found  that  the  unfolding  technique  under 
consideration  resulted  in  exposure  doses  which  were  substantial  overesti- 
mations  of  those  obtained  using  the  dosimeter.   This  overestimation 
appeared  to  be  systematic  and  it  was  concluded  that  the  most  probable 
source  of  this  error  was  in  the  normalization  of  the  response  functions 
and/or  in  the  generation  of  the  response  matrix. 


